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Abstract 

We expand a set of notions recently introduced providing the general setting 
for a universal representation of the quantum structure on which quantum infor- 
mation stands. The dynamical evolution process associated with generic quantum 
information manipulation is based on the (re)coupling theory of SU (2) angular 
momenta. Such scheme automatically incorporates all the essential features that 
make quantum information encoding much more efficient than classical: it is fully 
discrete; it deals with inherently entangled states, naturally endowed with a tensor 
product structure; it allows for generic encoding patterns. The model proposed 
can be thought of as the non-Boolean generalization of the quantum circuit model, 
with unitary gates expressed in terms of 3nj coefficients connecting inequivalent 
binary coupling schemes of n + 1 angular momentum variables, as well as Wigner 
rotations in the eigenspace of the total angular momentum. A crucial role is played 
by elementary j-gates {6j symbols) which satisfy algebraic identities that make 
the structure of the model similar to "state sum models " employed in discretizing 
Topological Quantum Field Theories and quantum gravity. The spin network sim- 
ulator can thus be viewed also as a Combinatorial QFT model for computation. 
The semiclassical limit (large j) is discussed. 

PACS numbers: 03.67.Lx, 03.65.Fd, ll.lO.Kk 

Keywords: Quantum Computation, Spin Networks, Topological Quantum 
Field Theory. 
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1 Introduction 



It is by-now a generally accepted fact that the laws of quantum theory pro- 
vide in principle a radically novel, and more powerful way to process informa- 
tion than any classically operated device p. In the past few years a big deal 
of activity has been devoted to devise and to implement schemes for taking 
actual advantage from such quantum extra power. In particular in Quantum 
Computation (QC) the states of a quantum system S are used for encoding 
information in such a way that the final state, obtained by the appropriate 
unitary time evolution of S, encodes the solution to a given computational 
problem. A system S with state-space Ti. (the Quantum Computer) supports 
universal QC if any unitary transformation U G W(7i) can be approximated 
with arbitrarily high accuracy by a sequence (the network) of simple unitaries 
(the gates) that the experimenter is supposed to be able to implement. The 
case in which iS is a multi-partite system is the most relevant, as it allows 
for entanglement, a unique quantum feature that is generally believed to be 
one of the crucial elements from which quantum speed-up (polynomial or 
possibly exponential) is generated Pj. 

In the above picture of QC the realization of the quantum network is achieved 
at the physical level by turning on and off external fields coupled to S as well 
as local interactions among the subsystems of S. In other words the ex- 
perimenter "owns" a basic set of time-dependent Hamiltonians that she/he 
activates at will to perform the necessary sequences of quantum logic gates. 

At variance with such a standard dynamical view of QC, more recently 
several authors considered geometrical and topological approaches [3]-[7]. 
The peculiarity of these proposals is far reaching: over the manifold C of 
quantum codewords one can have a trivial Hamiltonian, for example if |^ = 0, 
yet obtain nevertheless a non-trivial quantum evolution due to the existence 
of an underlying geometrical/topological global structure. The quantum 
gates - or parts thereof - in this latter case are realized in terms of operations 
having a purely geometrical/topological nature. Besides being conceptually 
intriguing on their own, these schemes have some built-in fault -tolerant fea- 
tures. This latter attractive characteristic stems out of the fact that often 
certain topological as well as geometrical quantities are inherently stable 
against local perturbations. This in turn allows for quantum information 
processing inherently stable against special classes of computational errors. 

There has been a number of proposals suggesting general conceptual 
schemes of interpretation of quantum computation. Most of them are in- 
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deed based on topological notions, even though this is not always explicitly 
stated. Among these, anyonic quantum computation 0, fermionic quantum 
computation fSj, localised modular functor quantum field computation [3], 
holonomic quantum computation [H], JHl have mostly attracted attention. 
However, such models appear to be simply different realizations of a unique 
conceptual scheme that incorporates all of them as particular instances (once 
one focuses on their suitable "discretized " counterparts). 

We propose here, expanding a set of notions introduced in ^T], a general 
setting for a universal representation of the very quantum structure on which 
quantum information stands. The associated dynamical evolution process, 
giving rise to information manipulation, is based on the (re) coupling theory of 
SU{2) angular momenta (see [12], [ISl Topic 12, [E]). The scheme automat- 
ically incorporates all the essential features that make quantum information 
encoding so much more efficient than classical: it is fully discrete (both for 
space-like and time-like variables); it deals with inherently entangled states, 
naturally endowed with a (non associative) tensor product structure; it al- 
lows for generic encoding patterns. The minimal set of requirements listed by 
Feynman ^5] as essential for the proper characterisation of an efficient quan- 
tum simulator is automatically satisfied: i) locality of interactions; ii) number 
of "computer " elements proportional to a function which is at most polyno- 
mial in the space-time volume of the physical system; iii) time discreteness 
(time is itself "simulated" in the computer by computational steps). 

Key element of our argument is the fact that all such basic features are 
typical of spin networks. It should be enphasized that by spin networks we 
mean here - contrary to what happens in solid state physics, but somewhat in 
the spirit of combinatorial approach to quantum space-time representation 
[TB] - graphs the node and edge sets of which can be labelled by quantum 
numbers associated with SU{2) irreducible representations and by SU{2) 
recoupling coefficients, respectively. For this reason spin networks can be 
thought of as an ideal candidate conceptual framework for dealing with ten- 
sorial transformations and topological effects in groups of observables. The 
idea is to exploit to their full extent the discreteness hypotheses ii) and iii), 
by modelling the computational space in terms of a set of combinatorial and 
topological rules that mimic space-time features in a way that automatically 
includes quantum mechanics. 

The model proposed can be thought of as a non-Boolean generalization 
of the quantum circuit model, with unitary gates expressed in terms of: a) 
recoupling coefficients (3nj symbols) between inequivalent binary coupling 
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schemes of = n + 1 S'?7(2)-angular momentum variables (j-gates); b) 
Wigner rotations in the eigenspace of the total angular momentum (M- 
gates). These basic ingredients of the spin network simulator, namely com- 
putational Hilbert spaces and admissible elementary gates, are discussed in 
details in Sections 2 and 3, respectively. The picture does contain the Boolean 
case as the particular case when all angular momenta are spin i. 

In Section 4 both the architecture and the computational capabilities of 
the simulator are described in full extent. On the kinematical side (Section 
4.1), the computational space is shown to be modeled as an S'[/(2)-decorated 
graph or, more precisely, as a fiber space structure over a discrete base space - 
the Rotation graph - which encodes all possible computational Hilbert spaces 
as well as gates for any fixed number of incoming angular momenta. A 
crucial role is played by elementary j-gates (Racah transforms, related to 6j 
symbols of SU{2)) which satisfy suitable algebraic identities [14] and make 
the structure of the model similar to "state sum models" employed in dis- 
cretizing Topological Quantum Field Theories (TQFT) and quantum gravity 
(see Section 5). In Section 4.2, after a discussion of the hypotheses of Feyn- 
man ^|, general circuit-type computation processes on the spin network are 
described and classified into computing classes. Virtual, polylocal Hamilto- 
nians are generated by the simulator which evolves in an intrinsic discrete 
time variable. Section 4.3 deals with questions in (quantum) computational 
complexity which turn out to be closely related to graph combinatorics. We 
argue that our new framework, when implemented on the basis of explicit 
encoding schemes |l7j, could be suitable to handle "combinatorially hard " 
problems more efficiently than any classical machine. 

The key ingredient of Section 5 is the Ponzano-Regge asymptotic formula 
[TH] for the 6j symbol which plays a twofold role. On the one hand, it gives 
a precise meaning to the semiclassical analog of a j-gate providing, together 
with the asymptotics of Wigner rotation matrices |13j (Topic 9), the notion 
of "approaching a classical simulator" out of a quantum one. Moreover, as 
stated in jTTj, the whole conceptual scheme of spin network computing can 
be reformulated in terms of density matrix formalism, namely resorting no 
longer to sharp eigenstates of angular momenta but rather to generalised 
multipole moments (see e.g. [19j, Ch. 7.7). Such generahsations can be 
summarised in the following diagram 
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SEMICLASSICAL SIMULATOR 



'll (large quantum numbers) 



SPIN NETWORK QU-SIMULATOR 



J| (extension to) -J| 
MIXED STATES COMPUTING MACHINE 



On the other hand, the Ponzano-Regge asymptotic formula opens the in- 
triguing possibility of bridging the quantum theory of angular momenta to 
Euclidean gravity in dimension three. More precisely, a state sum functional 
for triangulated 3-dimensional space-time manifolds, built up by associating 
a 6j symbol with each tetrahedron, is shown to correspond, in the asymptotic 
limit, to the semiclassical partition function of gravity ^H] with a classical 
action representing the discretized counterpart of the Einstein-Hilbert ac- 
tion of general relativity |5(T. Since quantum gravity in dimension three is 
strictly related to a TQFT with an SU (2) Chern-Simons-type action (see 
references quoted at the end of Section 5 and at the beginning of Section 
6) we present in this part of the paper some known results concerning spin 
networks viewed as "Combinatorial Quantum Field Theories", to be inter- 
preted as discretized versions of TQFTs based on S'?7(2)-decorated triangu- 
lations. At the end of Section 5 we compare the discrete partition functions 
of Ponzano-Regge gravity with the partition functions introduced in Section 
4.2 in connection with the simulator's dynamics. We conclude that a spin 
network simulator working by switching on j-gates acting on states with 
incoming spins is able to simulate some subclasses of triangulated surfaces 
and possibly some subclasses of triangulated 3-manifolds but, since partition 
functions for quantum field theories must be "sums over all configurations" 
(apart from regularization) , we cannot infer the possibility of fully simulat- 
ing Combinatorial QFT (unless we take some sort of termodynamical limit 
for N ^ oo which does not sound good when dealing with quantum circuit 
schemes for computation). 

We start Section 6 (Spin network and topological quantum computa- 
tion) by reviewing some basic definitions on TQFTs. Section 6.1 addresses 
holonomic quantum computation and we give indications that the discrete 
setting developed in Section 4 could support also such kind of computational 
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processes. In Section 6.2 we compare the spin network approach with the 
approach of Freedman and collaborators We provide a (not unique) 
mapping between the spin network and the modular functor approach by in- 
troducing combinatorial marked 2-disks which display localised interactions 
between spins. The algebraic structures of the two approaches, summarised 
in the Yang-Baxter identity for the standard topological one and in the 
(hexagon + pentagon) identities for the spin network, suggest that the parti- 
tion functions of the two models are related to each other in the same way as 
the regularized version of Ponzano-Regge functional corresponds to a double 
Chern-Simons partition function. 

In Appendix A we present the graph-theoretical rationale underlying spin 
network combinatorics by collecting results spread over a number of refer- 
ences in discrete mathematics, binary couplings and recoupling theory of 
angular momenta, complexity theory. Binary coupling trees are defined in 
Appendix Al, Twist-Rotation and Rotation graphs in A2 and some results 
in (classical) combinatorial complexity theory are summarised in A3. Ap- 
pendix B contains (standard) technical results concerning the composition 
of Wigner rotation matrices (Bl) and U-rotation matrices (B2) which are 
employed in Section 3.2 (M-gates). 

The twofold possible interpretation of the deliberately ambiguous title we 
have chosen for the paper should have become clear at this point: 

- on the one hand, spin networks are computing devices supporting simula- 
tions of the dynamical behaviour of composite quantum systems described 
in terms of pure angular momentum eigenstates; 

- such computing devices, on the other hand, are able to simulate classes of 
extended geometrical objects modeled as spin networks. 

We may summarise the content of the paper in the following diagram, 
where the spin network simulator may be viewed both as a generalised quan- 
tum circuit and as a Combinatorial QFT model for computation. The stan- 
dard Boolean quantum circuit (shown to be equivalent to the topological 
approach [T) is a particular case of this general scheme for computation. To 
complete the picture, the combinatorial approach can be suitably mapped 
into the purely topological one as discussed in Section 6.2. 

We plan to develop in the next future the upper connection which points 
toward "quantum automata" since our framework seems quite promising to 
address such issues like quantum languages and grammars, quantum encod- 
ing ^7] and quantum complexity classes of algorithms, naturally related here 
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to enumerative combinatorics of graphs. 



BOOLEAN CIRCUIT 



TOPOLOGICAL QFT 



QU-AUTOMATA 



GENERALISED 
QU-CIRCUIT 



SPIN NETWORK 
QU-SIMULATOR 



COMBINATORIAL 
QFT 



2 Computational Hilbert spaces 

Following (Topic 12) let us consider N = n + 1 mutually commuting 
angular momentum operators of the algebra of SU (2) 

Jl; J2; J3)---)Jn+l = {J«} 

and the corresponding components 

{Ji{z)}i=l,2,...,n+1 

along the quantization axis. For each i = 1, 2, . . . ,n + 1 the simultaneous 
eigenstates of the complete sets Jj and Ji(z) are 
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J*^ = ji {ji + 1) Ijiirii) 

Ji(z) \jimi) = TTli liiUli), (1) 

where we set ^ = 1 and the eigenvalues range over 

, _ 1 3 

—ji < "^j < ji (integer steps). (2) 

Denoting by 

= span { Ijirrii) } 

the {2ji + 1) -dimensional Hilbert space supporting the jj-th irreducible rep- 
resentation of SU{2), the tensor product 

= span { I ji nil) ® • • • ® | j„+i } (3) 

represents the simultaneous eigenspace of the 2(n + 1) operators {3j ; Ji(z)} 
and may be used e.g. to describe the state of = {n + 1) kinematically 
independent particles. By setting ji = j2 = • • • = jn+i = 1/2 in (jH)) we would 
get ®^ C^, namely the iV-qubits space of the Boolean quantum circuit model. 

To address the interacting case, we make explicit the basic assumptions 
on which the model of simulator discussed in this paper relies: 

1. The simulator computational states are suitable pure (n + l)-angular 
momenta states (to be described below) in Wigner-coupled Hilbert 
spaces of the total angular momentum operator 

Jl + Js + J3 + . . . + Jn+l = J (4) 

and of its projection J^. The corresponding quantum numbers, J and 
M = mi + 1712 + ... + rnn+i {—J < M < J in integer steps), la- 
bel the resulting representation spaces (which are referred to as JM- 
representations). 
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2. Interactions are not fixed a priori but rather are generated by the 
simulator itself as polylocal virtual Hamiltonian operators (see Section 



3. The admissible interactions are modeled on (a finite number of) com- 
binations of the following basic types: 

• binary couplings of the computational states which involve only 
the spin quantum numbers; 

• actions of rotation operators over states in JM-representations 
involving the total magnetic quantum number and depending on 
continuous sets of parameters {e.g. Euler angles). 

Given J as in (@)), the simultaneous eigenspace of the operators and 
(namely a JM-representation space for any fixed J) 



turns out to be degenerate. The degeneracy is partially removed by noticing 
that Jf, 3l, . . ., Jn^i commute with and Jz and thus ji, j2, ■■■,jn+i 
still are good quantum numbers (while the individual mi,m2, . . . ,mn+i are 
not). The ket vectors \J M) in (0) can be rewritten for the moment as 
lii; J2, • • • ,jn+i', J M), namely in terms of (n + 1) + 2 quantum numbers. 

The complete removal of the degeneracy would be achieved by introducing 
a new set of {n — 1) Hermitean operators - commuting with each other and 
with the previous ones - in order to get a total amount of 2(n + 1) quantum 
numbers (this number equals the number of operators needed to specify the 
eigenstates in the factorized Hilbert space 0). The most effective way to 
reach the goal is to consider equivalently: 

• binary coupling schemes in the sequence 



4.2). 



J2|JM) = J(J + 1) \JM) 
Jz \JM) = M \ JM) 



(5) 



Jl + 3 2 + J3 + ••• + Jn+l 



J 



(6) 



• binary bracketings on the factorized Hilbert space 



® ® eg) 



(7) 
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As a simple example, consider the case (n + 1) = 3: the binary coupling 
schemes are (Ji + J2) + J3 = J; Ji + (J2 + J3) = J; (Ji + J3) + J2 = J, 
and the corresponding binary bracketings are ® {W'^ ® 

W^)] (n^' ® W'') ® respectively. 

The crucial point is to realize that each binary coupling in ® generates 
- by using the Clebsch-Gordan series of SU{2) - an intermediate angu- 
lar momentum operator whose quantum number will be added to the set 
{ji, j2, . . . ,jn+i', J M}. In the case (n + 1) =3 the first coupling scheme 
(Ji + J2) + J3 = J splits into 

(Ji + J2) = J12 

J12 + J3 = J (8) 

with \ji — j2\ < ji2 < ji + 32 and J = ju + js, while the second coupling 
scheme Ji + (J2 + J3) = J splits into 

(J2 + J3) = J23 

Ji + J23 = J (9) 
with \j2 -jal < j23 < 32 + h and J = ji + jss- 

The reformulation of the same example in terms of binary bracketings on ((Tj) 
leads to the expressions 

{W^®{W^®W')i,,)j (10) 

where the inner brackets have been labeled by the q-numbers associated with 
the corresponding intermediate angular momentum operators and we have 
added an overall bracket labeled by the total J. Note that these Hilbert 
spaces, although isomorphic, are not identical since they actually corre- 
spond to (partially) different complete sets of physical observables, namely 
{Jj, J|, J^2; J35 J^) Jz] and {J^, J2, J3, JI3, J^, Jz] respectively (in partic- 
ular, J^2 and J23 cannot be measured simultaneously). On the mathematical 
side this remark reflects the fact that the tensor product is not an asso- 
ciative operation. 
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Coming to the general counting argument explained in Appendix 

Al (based on combinatorics of rooted labeled binary trees) shows that the 
number of binary bracketings one can accommodate on the (?T,+l)-fold tensor 
product in ((7j) is (n — 1), plus the external bracket (. . .)j. Thus we get exactly 
the number of intermediate angular momenta operators we need to remove 
completely the degeneracy in the JM space (0). More precisely, given a 
particular binary bracketing structure (for the moment we can think of an 
ordered sequence of incoming angular momenta {Ji, J2, J3, • • ., Jn+i}), we 
get a unique set of (ordered) mutually commuting operators denoted by 

Ki, K2, K3, ...,K„_i (11) 

with quantum numbers fci, ^2? k^i ■ ■ ■ ikn-ii each running over a suitable 
finite range (c/r. (jSI) and Q). 

An explicit example of one possible bracketing structure is given by 

(. . . ® ® n^%, ® . . . ® n^")k^_, ® w-^')k^_,)j (12) 

where the incoming angular momenta are coupled sequentially. To denote 
the basis vectors belonging to such a space we could write either 

J2, j3,---,Jn+i;fci,...,fc„_i; JM), -J <M<J} (13) 
(where both the sequences of quantum numbers j's and fc's are ordered) or 

{ I (. . . ((((ji, 32)k, ,j3)k„ ■ ■ .), Jn+i)fc„_Jj ■,JM),-J<M< J}, (14) 

where the binary bracketing structure (. ® •)k of (fT^ has been exactly 
transferred on the string of quantum numbers inside the symbol | . . .) as 

( ■ ; Ofc • 

However, when dealing with other types of binary bracketing structures 
none of the above notations turns out to be well suited. First, we would 
like to consider any permutation of the incoming angular momenta quantum 
numbers {ji, j2, js, • • • ,jn+i} and not just a sequence with a fixed ordering. 
It is clear that the notation adopted in (jl3|) is not fiexible in this respect. 
Secondly, we have to get rid of all possible binary arrangements of the in- 
coming variables and at the same time of the resulting (partially ordered) set 
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of intermediate fc's which appears in the example (jl4|) as subscripts under 
brackets. 

According to the above remarks we shall denote from now on a binary 
coupled basis of (n + 1) angular momenta in the JM-representation (and 
the corresponding Hilbert space) as 

{ I [Jl, J2, J3, • • • , Jn+l]' ; klkl..., kl_, ■ JM), -J<M< J} 

= Hi (b) = span{ \ b-JM)n} , (15) 

where the string inside [ji, j2, js, ■ ■ ■ , jn+i]^ is not necessarily an ordered one, 
b indicates the current binary bracketing structure and the fc's are uniquely 
associated with the chain of pairwise couplings given by b. 

As explained in details in Appendix Al, the combinatorial structure un- 
derlying the above computational Hilbert spaces (for fixed n and for any J) 
is provided by rooted labeled binary trees with: 

• (n + 1) terminal nodes (or leaves) labeled by {ji, j2, js, • • • ,jn+i}] 

• (n — 1) internal nodes labeled by {ki, k2, fcs, . . . , kn-i}; 

• the root, labeled by J. 

For instance, the binary coupling tree structure corresponding to (jl2j) and 
flT^ is depicted in Fig. 12 of Appendix Al. 

As shown in the following Section, processing information through changes 
of binary coupling schemes of the Hilbert spaces ()15|) involves only spin quan- 
tum numbers (c/r. the first assumption in point 3.) and can be modeled on 
other types of combinatorial structures, Twist-Rotation and Rotation graphs 
as discussed in Section 4.1. 

On the other hand, the quantum number M can be activated by noticing 
that a vector in TC^ (b) is expressed in terms of the basis (fT3j) as 

J 

\b;J)n = Yl ^'^(^i\b;JM)n, (16) 

M=-J 

where ^^^af^ (M = — J, — J + 1,...,J — 1,J) represent the (complex) com- 
ponents of the vector. The natural unitary transformations acting on the 
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M-dependence of basis vectors and components (Wigner rotation operators 
^'m M'i^^Pl) introduced in Section 3.2) are not discrete, but rather depend 
on continuous parameters as claimed in the second assumption of 3.. 

Remark 2.1. A more fundamental question arising in connection with 
assumptions 1., 2., and 3. is the following: we could have included the 
factorized states (jH)) too in the set of computational states of the simulator, 
e.g. as non interacting input states. Then any binary coupled state (fT^ could 
be generated by means of a suitable sequence of Clebsch-Gordan coefficients 
representing a unitary transformation depending on the whole set of quantum 
numbers {ji, j2, ■ ■ ■ , jn+i, mi, m2, . . . , m„+i; ki, k2, . . . , kn-i, J, M} (see [H], 
Topic 12). For instance, a vector in the binary basis associated with (jSJ is 
obtained from the factorized basis through 

I((ji,j2),i,, j3)j; JM) 

= E C/^..^.m,CjX}:m, Ijimi)®|j2m2)®|j3m3), (17) 

mi m2 ma 

where there appear two Clebsch-Gordan coefficients. This expression can 
obviously be inverted, providing us with one factorized basis vector in terms 
of a combination of binary coupled states in the JM-representation. 

However, such kind of procedure - namely going through the factorized 
basis at each step of "computation " - is not satisfactory in many respects 
and certainly very cumbersome when dealing with many spins. Even if we 
could assume that the input state of the simulator is indeed a factorized, 
non-interacting one (and thus a transformation of the type (fTTj) is certainly 
allowed), any further unitary computing step would force us to deal with 
(pure angular momenta) coupled states. Now, quantum theory of angular 
momentum gives us a powerful framework to manipulate the binary coupled 
states of the type we are considering: other kinds of couplings {e.g. sym- 
metrical) cannot be efficiently manipulated since the degeneracy in the total 
JM-representation space cannot be completely removed for N = n + 1 > 3 
(see the remark at the end of [iSj, Topic 12). 

Summing up, a model of simulator whose computational states are binary 
coupled ones actually relies on the most exhaustive theory we have at our dis- 
posal in order to deal with many-angular momenta quantum systems whose 
interactions are modeled (or can be well approximated) by polylocal two- 
body-type interactions. It is worth noting that two-body interactions were 
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shown to be sufficient to implement universal quantum computation in the 
decoherence-free subspaces approach according to j2I] (see other references 
therein). A 



3 Gates 
3.1 j— gates 

By j-gates we mean unitary transformations on the computational Hilbert 
spaces flT^ which act on the set of the spin variables {ji, j2, ja, • • • , jn+i, 
fci, ^2, ^3, • • • , kn-i} of the eigenstates without changing the quantum num- 
bers J and M. 

According to the recoupling theory of angular momenta ^21; (Topic 
12), (see also [22]; 123) the most general unitary transformation between 
two computational states characterised by different binary coupling schemes 
b and b' 

I [jl, i2, is, • • • , Jn+l] '-^klkl..., kl_^ ■,JM)^ 

^ I [Ji, J2, J3, • • • , Jn+l]"' ; fcf , . . . , ; JM ) (18) 
is a recoupling coefficient of SU{2) (or 3nj-symbol) denoted by 



U 



3nj 



'^1 ■ ■ ■ 

'^1 ■ ■ ■ li'r, 



Wan, [b;b'], (19) 

where the variables {j}, {k'^}, {k^'}, J-, M appearing in states (fTHj) have 
been partially or totally dropped (c/r. the shorthand notation for Hilbert 
spaces and state vectors introduced in (US))). According to this definition, 
\Usnj [b; b']p represents the probability that a quantum system prepared in 
the state |b; JM)n is measured in the state |b'; JM)^. 

Note that (fT^ is a reduced tensor operator, namely the magnetic quan- 
tum numbers do not appear at all in its expression in view of the Wigner- 
Eckart theorem ^H]- More explicitly, by looking at the recoupling between 
two states of the type (fTH|) -but with different M, M'- we would get 



( [ji, . . . , j.+i]"' ; fcf , . . . , kl_, ■ JM' I [ji, . . . , jn+i]"; . . . , /cti ; JM) 
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= 6MM'{[...f;kl',...,k'^_,; JM'||[...]^fc^...,fc^l; JM), (20) 

where S is the Kronecker delta and the symbol (. . . 1 1 . . .) denotes the reduced 
operator (which coincides with U^nj [b; i>'] in ()19|1 for M = M'). 

For each fixed n there exist inequivalent types of 3nj symbols (up to de- 
generate cases and phase transforms, see below): one 6j symbol, one 9j, two 
12j symbols, five 15j's, eighteen 18j's etc. ^2)- The recoupling coefficients 
classified as type I and type II may be expressed through single sums of prod- 
ucts of 6j symbols, while types III, IV, V, . . . (appearing for n > 5) may 
be represented either by single sums of more complex products of 6j and 9j 
symbols or by multiple sums of products of symbols of lower orders. In the 
present context we do not really need any such complicated expressions since 
we take advantage of the results illustrated in details in (Topic 12) and 
collected in the following: 

Theorem (Biedenharn-Louck). 
For each n any U^nj \^'-, is the composition of (a finite number of) two 
elementary unitary transformations, namely 

• Racah transform 

n : \...{{ah)dC)f...-JM) ^ \...{a{bc)e)f...;JM), (21) 

• Phase transform 

$ : \ ...{ab)c...;JM) ^ \ . . . {ba)c . . . ; JM), (22) 

where we are using Latin letters a,b,c, . . . to denote both incoming (j's in 
the previous notation) and intermediate (k's) spin quantum numbers. 
The explicit expression of (jTIj) reads 

|(a(6c)e)/;M) = 
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where there appears the Racah-Wigner 6j symbol of SU{2) and the weights 
{2d + 1), (2e + 1) are the dimensions of the irreps labeled by d and e, respec- 
tively. 

Recall that the 6j symbol may be expressed in turn as a sum over magnetic 
quantum numbers of the product of four Clebsch-Gordan coefficients with 
entries in the set {a,b,c,d,e, f; ma,mh,mc,md,me,mf}, where (—a < 
nia < a in integer steps) is the magnetic quantum number associated with 
the spin variable a (and similarly for the others) |14j. The numerical value 
of the 6j symbol depends on normalization: we are tacitly assuming through 
the whole paper the standard Condon-Shortley conventions. 
Finally, the phase transform (|22j) reads 

\...{ab),...-JM) = {-ir+''~'\...{ba),...-JM). A (24) 

Summing up, the unitary transformations {U^nj [b; b'] } connecting any 
pair of computational states |b JM)n, |b; JM)n are identified as j-gates in 
the present computational context and the theorem stated above naturally 
single out two types of universal elementary j-gaXes. In particular, the Racah 
transform switches the current coupling scheme to a physical inequivalent one 
(c/r. the explicit expression (j22I))- 

As a further remark note that the elementary unitary j-gates ()21|) and 
fl22|) are in one-to-one correspondence with elementary topological opera- 
tions between pairs of rooted labeled binary trees. Rotations and Twists 
respectively (c/r. Appendix A2 and pay attention to the fact that "rotation" 
will be used in the following part of the present Section in a totally different 
way). 

3.2 M-gates 

As anticipated in Section 1, M-gates are implemented -inside each computa- 
tional Hilbert space H'^ (b) introduced in (fT^ - by the Wigner D-functions 
DIj j^p, namely the matrix elements of the (unitary) rotation operator T)^^, 
in the JM representation. By using Euler angles a, (3, 7 to specify the 
rotation, the Wigner Z)-functions -or W-rotation matrices- are defined by 

( JM («/37) I J'M') = bj,r («/?7), (25) 

where < a < 27r; 0</3<7ror27r</3< Svr; < 7 < 27r and we adopt 
in this Section the notations of |14j (see also |19j for a general treatment of 
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SU{2) representation theory). A W-matrix can be always expressed as 

D'mm' («/57) = e-^^" di^, e-^^'\ (26) 

where c^mm' iP) the reduced W -matrix. 

Generally speaking, W-rotation matrices appears when we look at trans- 
formations of the wave function of a quantum mechanical system under a 
coordinate rotation. Given a basis vector of (b) as in (P3|) . its depen- 
dence on polar coordinates {9, (p) is better described if we change our previous 
notation according to 

|b;JM)„ ^ ("Vifl^,*/)) (27) 

where we dropped the subscript n for simplicity. Then the action of a W- 
matrix on any such basis vector is expressed by 

+J 

("^^i,, (^',0') = E ^''Vm(^,0) ^mm'Mt), (28) 

M=~J 

where {9, 0) and (6'', 0') are polar angles in the initial and rotated coordinate 
systems, respectively. The action of a rotation on the components of vectors 
(c/r. the decomposition written in involves the Hermitean conjugate 

D ppj^j of the corresponding Df^j^^,, namely ^''^ a'/j, = ^D j^^j^j ^''^ afj, where 
the sum is over M as in (|28p . 

It is worth to recall that in the general case (A^ any integer and ji,j2,..., 
jn+i chosen in {0, 1/2, 1,3/2, . . .}) the reducible (2J+l)x (2J+1) W-rotation 
matrix Dj^jj^p will admit a block diagonal decomposition into irreducible 
rotation matrices of lower ranks. On the other hand, by changing N and/or 
the values of the incoming spins, the resulting -D]^^, will decompose into 
different elementary blocks. From the computational viewpoint this provides 
a more general notion of universal set of elementary M-gates than that 
currently adopted in (Boolean) quantum information schemes, typically given 
in terms of 2 x 2 and/or 4x4 unitary matrices [21]. 

However, as described in details in Appendix Bl, each matrix element 
of any D'lij^, can be factorized in a well-defined way, and this procedure 
is independent of the binary bracketing structure of (b). The explicit 
expression of such factorisation given in ()135|) can be written symbolically as 
in (jl36|) . namely 
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N 

E n (C^o, D'H^Pl) C:i,,^), (29) 

{m,m'} i=l 

where we dropped the matrix indices M, M' ,mi, m[ on the W-matrices D"^ , 
D^^ {i = 1,2, ... ,N = n + 1) and similarly C*,'s are Clebsch-Gordan co- 
efficients with m-type entries omitted. The summation is over all magnetic 
quantum numbers of the angular momentum operators {Jj}, while {kj} are 
spin quantum numbers associated with the intermediate operators defined in 
(USUI)- 

In ()29|) there appear N {= ^ of incoming spins) factors, each containing 
a W-matrix in the irreducible j^-th representation of dimension (2jj + 1), 
and a total amount of 2N Clebsch-Gordan coefficients. As explained at the 
end of Appendix Bl, each matrix element of D^^ may be further factorized 
into the (sum of) product of 2jj W-matrices in the fundamental j = 1/2 

representation of SU{2). The explicit form of D^^, {m,m' G {\, —\}) 
reads 



1 /e— /2 cos(/3/2) -e-^-lHm{P /2) e^^l^\ 

^mm' l"P7J - gia/2gin(^/2) e^"/^ cos(/5/2) 6^^/^ )■ ^"^^^ 

As a consequence of the above remarks we conclude that the "elementary 
factors" appearing in the right hand side of needed to determine one 
matrix element of D'^ , namely D^j^, (a/37) fo^^ some MM', are 



2N C-G coefficients 

N 

2J= 2J2 Ji W-matrices D^^, (31) 
1=1 

and the number of factors one needs in order to evaluate the whole D'^ 
amounts to [2( J + A^)]'^^'^+^) . These estimates represent in fact upper bounds 
on the number of factors, since we may reduce the number of elementary W- 
matrices by employing some 3x3 matrices of the j = 1 irrep. Moreover, 
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by considering a purely fermionic (bosonic) symmetric A^-multiplet, the ex- 
pression (|29p does not contain C-G coefficients anymore and the number of 
elementary factors to be taken into account in ()31|) is simply 2J = N [J = N, 
respectively), as discussed in Appendix Bl. 

Remark 3.1. Our framework turns out to be much richer than the 
Boolean case, but contains the usual Boolean gates as particular examples. 

Note also that D^^, is an elementary gate in any situation; it is also universal 
for the two particular cases discussed above. In particular, the N spin-i case 
is compatible with the scheme proposed in j2Tj, Sect. VII. 

Remark 3.2. In each representation labeled by J the W-matrices form 
a group under multiplication, namely 

(«i/?i7i) i:>^ («2/?272) =D'{a(3^), (32) 

where (a/?7) are related to (ai/?i7i) and (0:2/5272) by quite involved expres- 
sions (see [14J Ch. 1.4.7). 

Remark 3.3. Instead of W-matrices we could have used other parametri- 
sations for rotations, for instance the matrices t^/f^v/' defined in terms of ro- 
tation axis and rotation angle. For completeness we collect in Appendix B2 
some standard formulas relating these two types of transformations together 
with some explicit examples ^1]. ▲ 



4 Spin Network Quantum Circuit 

By exploiting the basic ingredients introduced in the previous sections (com- 
putational Hilbert spaces, j-gates and M-gates) we present here the struc- 
tural setting of a quantum simulator DJl -the spin network simulator- mod- 
eled as a generalised {i.e. not Boolean) quantum circuit model. In a broader 
sense such a computing machine could be reinterpreted as a concrete realiza- 
tion of what should be a Quantum Automaton (see e.g. [23, ESI), namely a 
theoretical framework able to deal consistently with quantum languages and 
grammars [T7j . 
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4.1 Combinatorial kinematics 

The computational space of the simulator DJl turns out to be modeled on 
an S'f/(2)-decorated graph (or spin network, according to Penrose's similar 
structures introduced in [TE]). 

For each fixed n the underlying network structure is denoted for the 
moment by Gn{V,E), where V and E are the vertex and edge sets of 
the graph, respectively. Both the vertices and the edges (arcs connect- 
ing pairs of vertices) of Gn(y,E) are decorated by algebraic objects from 
S'f/(2)-representation theory we introduced previously. The vertices are 
in one-to-one correspondence with the set of computational Hilbert spaces 
^n(ii)i2, • • • , jn+i; h,---, K-i, JM) = Hi{ b) introduced in ^ 

vertex set 1/ < — > (33) 

and since each 'Hnih) has dimension (2J + 1) over C there exists one isomor- 
phism 

Hi{b) =1, C^'^+i (34) 

for each admissible binary coupling scheme b of {n+1) incoming spins. Thus 
each vertex in the spin network is decorated with a copy of C^"^"*"^ 

C^-^+i ^ v{b) G V, (35) 

where each v & V is labeled by the unique corresponding b. 

The construction of the graph proceeds by establishing the connections 
between vertices. The edge set E = {e} of Gn(y,E) is a subset of the 
Cartesian product (V x V) selected by the action of elementary j-gates. 
More precisely, an (undirected) arc between two vertices f (b) and v{b') 

e(b,b') = {v{b),v{b')) G (VxV) (36) 

exists if, and only if, the underlying Hilbert spaces are related to each other by 
one of elementary unitary operations defined in Biedenharn-Louck Theorem 
of Section 3.1. Note that at the kinematical level the resulting decorated 
edges in the net are to be considered undirected since any such operation is 
invertible. 

This general combinatorial structure can be concretized in two different 
ways, depending on what types of binary bracketing schemes -associated 
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with the vertex set of G„(V, E) as stated in (jH^ - we are going to consider as 
inequivalent (we refer to Appendix A2 for more details on the combinatorics 
of the resuhing graphs). If we distinguish pairs b, b' which differ either by a 
Racah transform TZ ()21|) or by a Phase transform $ ()22|) . then the structure 
we get is a Twist -Rotation graph 

G^^(F,E) = &n{V,E) (37) 

with 

&n D V = {v{b), all admissible b} 

^ E = {e(b, b') ^ 7^, $}. (38) 

For each n the Twist-Rotation graph is a regular, cubic (trivalent) graph rep- 
resenting pictorially all possible types of 3nj symbols introduced in Section 
3.1 and extensively analysed in literature (c/r. [221, IS], |2Z|)- In Fig- 
21 of Appendix A2 the Twist-Rotation graph ©3 (binary coupling schemes 
of four angular momenta and associated 9j symbols) is depicted. 

Since quantum states which differ by phase transformations give rise to 
the same physical probabilities of observables, it seems quite natural to con- 
sider a combinatorial structure in which b, b' correspond to different vertices 
only if the states are related by a Racah transform (recall that in this case 
there is a change of one intermediate spin variable). This reduction opera- 
tion on the Twist-Rotation graph ()37p involves an equivalence relation to be 
imposed on the vertex set. Moreover, each of the surviving edges turns out 
to be associated with a Racah transform TZ, but there could appear some 
additional weights and/or phase factors in its explicit expression (j23I)- With 
these premises we define the Rotation graph 

G^iV,E) = eAV,E) (39) 

with 

0„ D V = {equivalence classes off (b)} 

(5n E = {e(b, b') TZ, mod weights/phases}. (40) 

Figure 1 
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A picture of ©3 is given in Fig. 1, while more details on the above construc- 
tion are collected in Appendix A2. Here we just recall that the Rotation 
graph lJ5„ is a regular (not planar) 2{n — l)-valent graph with a number of 
vertices given by 

\V\ = card\/ = (2n - 1)!! (41) 

where ()!! denotes the double factorial number. In the following we shall 
deal mainly with the Rotation graph structure, but we point out in advance 
that applications of such combinatorial machinery in the computing context 
work on the Twist-Rotation graph as well. 

On the spin network the action of a j-gate defined by a recoupling co- 
efficient Usnj [b] b'] introduced in (fTUI) is represented formally by a piecewise 
path in 0„ connecting the vertices v{b) and v{b') 

N '^snj [b;b'] 

v{b) > v{b). (42) 

Two state vectors [e.g. basis vectors, for simplicity, with the same value of 
the q-number M) in the Hilbert spaces (j^Hj) attached to the vertices labeled 
by b (the initial point) and b' (the terminal point) are denoted by 

|b; JM)rr = |in(b); JM)„ and |b'; JM)„ = |out(b'); JM)„ , (43) 

and can be interpreted consistently "in " and "out " as referring to an initial 
state (input) and a final state (output), respectively. 

To deal with the formal expression ()42|) we need some results from angular 
momenta recoupling theory summarised in the following 

Theorem (see ^2]; QH| Topic 12, and the original references therein). 

Consider all possible paths in (S^ connecting two states as in PHj). Then 
the transition probability amplitudes for any pair of paths with given end- 
points, say 7i and 72, are equal 

(in(b); JM|out(b'); JM)^, = (in(b); JM | out(b'); JM)^^ (44) 

since we can freely deform such paths one into each other. As a consequence 
of that, probabilities too turn out to be equal owing to the fact that 

\U3nj[b;b']\'^ = I (in(b); JM|out(b'); JM) |2 (45) 
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is actually the probability that a quantum system prepared in the state "in " 
will be measured in the state "out ". The proof of (j44j) relies on the existence 
of fundamental algebraic identities involving 6j symbols, namely 

• the Biedenharn-EUiott identity 

' ^ '\c d pj\e f qj\h a rj 
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• the Racah identity 

where the spin variables {a, 6, c, . . . , x} run over {0, |, 1, |, . . .} and must 
satisfy suitable triangular inequalities inside each 6j symbol (otherwise the 
symbol itself would vanish). The weight (2x + 1) is the dimension of the rep- 
resentation labeled by the quantum number x, the sum over x is constrained 
only by the triangular conditions quoted above and R in the phase factor of 
the first identity is the combination {a + h + c + d + e + f+p + q + r). Note 
that these identities, together with the orthogonality relation 

define uniquely the Racah-Wigner 6j symbol (considered as the hypergeo- 
metrical polynomial which generates the Askey hierarchy [2H])- 

Without entering into details about the proof of the above theorem, we 
collect below some remarks which should make the rationale of its proof as 
clear as possible. 

• The spin network 02 (see Fig. 22 in Appendix A2) is a closed loop 
of triangular shape: the vertices correspond to the three inequivalent 
binary couplings of (n + 1) = 3 spins and each edge is associated with 
a Racah transform. Then Racah identity (jTTj) ensures that we may 
implement the transition from one vertex ("in") to another ("out") 
traversing either the edge connecting them directly or the other two. 
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• The spin network 03 shown in Fig. 1 is characterised by triangular and 
pentagonal closed loops (bounding triangular and pentagonal plaquet- 
tes). Triangles are associated with Racah identity (jTrf) (as happens 
for (J52), while pentagonal plaquettes turns out to correspond to the 
Biedenharn-Elliott identity if the chosen path embraces two edges 
of a pentagon, one can freely deform it traversing the other three edges 
(and viceversa). By using both (j46|) and orthogonality relation (j48|) 
we can also deform a piece of path connecting two contiguous vertices 
into a path which touches the other four edges of the pentagon (and 
viceversa). Thus we conclude that the improvement of ()44|) for 03 
relies on all the three algebraic identities written above. 

• The spin networks 0„ (n > 3) display plaquettes with other types of 
polygonal boundaries. Each type of plaquette can be associated with 
a suitable algebraic identity which can be derived from the fundamen- 
tal ones by making use of the explicit expression for the 3nj symbol 
involved. The procedure for improving (jl^ goes on as in the previous 
cases. 

As a final comment on the equi-amplitudes of paths in 0„ under the action 
of pure j-gates, notice that such result holds as far as we are dealing with 
the computational space of the simulator at the kinematical level. When we 
shall ask the spin network to perform a computation by means of a sequence 
of j-gates such an invariance will be broken (as we are going to explain in 
the following section). ▲ 

In order to include in the combinatorial setting the action of M-gates 
(see Section 3.2) we may employ either the Wigner rotation matrices (I25|) 
parametrised by Euler angles or [/-matrices written in terms of rotation axis 
and rotation angle introduced in Appendix B2: here we agree to making use 
of the former ones as in the rest of the main text. By analogy with the formal 
expression (jl^ representing a j-gate on the spin network, we write formally 
the action of an M-gate on 0„ as 

vib) vib), (49) 

where we drop the matrix indices M, M' by using the operatorial notation 
as in ()25|). Since rotations do not alter the binary bracketing structure of 
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the computational Hilbert space 7i;( (b) we may activate an M-gate indepen- 
dently at any vertex v{b) C C5„,. To recover the explicit expression of the 
action of an M-gate we have to pick up a vector in 7i;((b), change back our 
notation as in (jTfj) (to make the angular dependence explicit) and finally re- 
cover the expression given in (j^Hj) . However, in order to have at our disposal 
a unified notation for states to be considered as "input " and "output " in a 
quantum circuit framework, we are forced here to use a hybrid notation by 
setting 

ni{b) = {span \b;e, 0; JM)„ }. (50) 

With this convention we write down the action of an M-gate (for a given 
choice of the parameters 7) on an input (basis) state as 

V'{a/3^) : |b;in(^, 0; M)J)„ ^ |b;out(^', 0'; M')J)„ (51) 
where the output state is 

J 

|b;out(5', 0'; M')J)„ = ^mm' Mt) 0; JM)„ (52) 

M=-J 

and the input basis state appears in the combination on the right-hand side 
with its particular M-label. 

The actions of both types of gates on the spin network can be visualised 
by looking at Fig. 2: we can move from one vertex to a different one along 
an edge as in (j42j) (without changing M ) or choose to perform a rotation 
()49|) inside the computational Hilbert space associated with a vertex. 

Figure 2 

According to (plj) and (fTBj) . a suitable unified notation for generic state vec- 
tors to be used in actual computations should be 

|b;^, 0;J)„ (53) 

and consequently input /output states in the particular cases (HHj) and (|^ 
have to be set in the form 

|in(b;^, 0); J)„; |out(b; 0'); J )„ (54) 
possibly with additional M-labels if basis vectors are considered. 
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Summing up, the kinematical ingredients of the spin network simulator 
based on the Rotation graph 0„(V, E) are 

E = {e(b, b')} elementary j-gates 

{V' : niib) ^ni{b)}^M-gates (55) 

where the discrete structure encoded in (V, E) is endowed with transforma- 
tions V"^ depending on both discrete and continuous parameters (c/r. the 
points discussed in Section 2). 

Remark 4.1. The combinatorial setting described above can be inter- 
preted as a fiber space structure {V,C^'^^\ SU{2y^) where 

• V = {f (b)} is the (discrete) base space; 

• C^'^^^ is the typical fiber, a copy of which is attached to each t>(b) 
through the isomorphisms given in ()34|): 

• SU{2y is the automorphism group of the fiber realized by the (2 J + 
1) X (2 J + 1) W-matrices which form a group under the composition 

law dig); 

• E = {e(b, b')} are arcs connecting pairs of contiguous vertices in the 
base space (c/r. but they may be also considered as mappings 

ivib),niib)) ^ iv{b'),niib')) (56) 

connecting each given decorated vertex to one of its nearest 2(n — 1) 
vertices. 

In such a vector bundle framework one could take advantage of the above 
transport prescriptions (W-matrices along the fiber and maps (j56|) along 
horizontal sections) to get a notion of "connection" in the total fiber space 
V X C2'^+^ as illustrated in more details in |29^. This remark opens the 
possibility of discussing relations between the spin network scheme and the 
holonomic q-computation approach. We shall come back on this point in 
Section 6.1. ▲ 
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4.2 Dynamics and computing 

The kinematical structure of the spin network 971 comphes with all the req- 
uisites of an universal q-simulator as defined by Feynman namely 

• locality, reflected in the binary bracketing structure of the computa- 
tional Hilbert spaces, which -together with the action of W-rotations 
- bears on the existence of local interactions; 

• discreteness of the computational space, reflected in the combinatorial 
structure of fi5„; 

• discreteness of time, to be discussed below; 

• universality, guaranteed by the properties of gates we described in Sec- 
tion 3: any unitary transformation operating on computational Hilbert 
spaces can be reconstructed by taking a finite sequence of Racah trans- 
forms (and possibly phases) followed by the application of a finite num- 
ber of W-rotations. 

Thus we have explicitly defined the class of "exact imitators" of any finite, 
discrete quantum system (described by pure angular momentum states) with 
no need of resorting to the notions coming from the (inherently classical) 
Boolean circuit theory. 

In order to describe the dynamical behavior of the spin network, we no- 
tice preliminarily that the rule to "move" from a state (say a vector (jSHj) ) 
to a nearest one have been already established: apply either one j-gate or 
one particular M-gate (for fixed a, /5, 7 ). Thus a natural discrete time unit, 
denoted by r, can be associated with one elementary step in such a cellu- 
lar automaton scheme. However, as pointed out by Feynman himself, this 
naive assumption can at best make the simulator to "imitate" time. If we 
pretend the spin network OK to "simulate " time, we have to go through a gen- 
uine space-time dynamics providing Hamiltonians and intrinsical evolution 
in actual time intervals. 

Generally speaking, the basic data to implement computation in a circuit 
model are an input state and a program V giving instructions to manipulate 
information stored in the machine states: output states must belong to the 
set of "accepted" states (if computation halts, as we are tacitally assuming). 
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In the spin network C5„ the choice of a particular program V is interpreted 
as the selection of a subset of unitary transformations 

V < — > {"dv} d {il} (57) 

among all the kinematically allowed {it}. Since we are going to deal with 
sequences of states we are forced to change again our last notation into 
a simplified one, namely 

|b; e, (j)- J)n ^ (58) 
Then a computation based on the program P, represented formally as 

|Oin)n ^ |Dout)n (59) 

is a collection of directed paths in C5„, all starting from the input state |t)in )n 
and ending in some accepted |t)out )n- By a directed path we mean a (time) 
ordered sequence of states 



|f in )n = |00 )n |0l )„ ^ ^ |D. )„ ^ ^ |0l )n = |Oout )n (60) 

where s = 0, 1, 2, . . . , L is the lexicographical labelling of the states along the 
given path and L is the length of the path, which turns out to be proportional 
to the time duration of the computation process L ■ t = T in units of the 
discrete time step r. The integer L characterising the particular directed 
path in (lUUj) represents the number of time-ordered elementary operations 
(computational steps) needed to get |t)out)n from |t)in)n by employing the 
program V. It should be clear that from the same input |t)in)n the program 
V could select different paths to get (possibly) different output states. For 
instance we may represent by 



|Oin )n=\%)n^K)n^ ' K)n ^ ' I^o'l' )n = jfout )n (61) 

another path of length L' in the collection ending in an accepted \^'ont )n- 
Each arrow in ()60|) or ()61|) stands for one of the elementary operations (Racah 
transforms, Wigner rotations) described in the previous sections (c/r. also 
Fig. 2). Using our current notation an elementary computational step is 
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either |t)s )„ ^ jo^+i )„ 

or|Os)„ > (62) 

and the expression of a particular unitary transformation M-p (of length L) 
in fl59|) turns out to be a well defined time-ordered composition 

ilp = Ul o Ul^i o ■ ■ ■ 0U2 o Wi , (63) 

where each W is given by one of the operations (jU^ . 

The framework developed so far is an extremely flexible and powerful 
circuit modelization of quantum computing and we shall examine later on in 
this section the range of different types of computations that can be actually 
carried out. However, to complete the dynamical setting based on the iden- 
tification made in ()57p and of a program V with a collection of directed 
paths in 0„, we have to call into play algorithms. Then the program V{A) 
to perform a particular algorithm A is the specification of a suitable directed 
path in the collection (jSHl), starting from a given |tiin)n and ending in one 
particular |Oout)n- The associated unitary transformation is denoted by 

V{A) ^ ^v(A) C {ii^} (64) 

and for the path itself we may use the same notation as in ()60|) . Alternatively, 
by broadening the meaning of the symbol V{A), we agree that it represents 
also an ordered sequence of labelings 

V{A) ^ {s = 0, 1, 2, . . . , L{V{A))} (65) 

which turns out to be in one-to-one correspondence with the states {|Os)n, 
s = 0,1,2, L{V{A))}. The time lapse required to get the output is 
nothing but the length of the path in units of the discrete time step r , 
namely 

T{P{A)) = L{V{A))-T. (66) 

A circuit-type computation process in 0„ based on a program V perform- 
ing the algorithm A (which could be formally written as in (|59|)) is actually 
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represented by the expectation value of the unitary operator (jMj) evaluated 
between the given input and output 

(Oout I ^V{A) I ^in )n- (67) 

This quantity gives the physical transition probability amplitude to get |Dout )n 
by acting with ilp(^) on |t)in and obviously its square modulus is the quan- 
tum probability to be assigned to the corresponding computation. By taking 
advantage of the possibility of decomposing itp(^) uniquely into an ordered 
sequence of elementary operators (gates), (jU7|) becomes 

L-1 

(Oout I *^V{A) I Oin )n = VW {^s+l \ i^s,s+l \ )n \ v{A) (68) 

with L = L(V{A)) for short. The symbol [ \ p(a) denotes the ordered prod- 
uct along the path V{A) and provides a sort of superselection rule which 
induces destructive interference of the forbidden {i.e., not leading to the cor- 
rect result) paths in Each elementary operation introduced in the generic 
expression (jUHj) is now better denoted by lAs^s+i to stress its "one-step " char- 
acter with respect to computation. Consequently, each elementary transfer 
matrix in (jUHjl turns out to be associated with a local Hamiltonian operator 
arising from 

(Os+i I W^.s+i I )n = exp{iiH„(s, s + 1) -r} (69) 

and representing the unitary evolution of the simulator in one unit of its 
intrinsic time variable (s = 0, 1, 2, . . . L{V{A))). We indicate with the short- 
hand notation (s, s + 1) the dependence of H„ on its variables to make clear 
the local nature of this operator with respect to the computational space (J5„ 
(a more detailed description should involve the quantum numbers character- 
ising both states and operation). When (jU^ is inserted in such virtual 
Hamiltonians generally do not commute with each other but nonetheless the 
whole computational process may be identified with a well defined unitary 
evolution of the simulator in the internal time interval given in . 

The above remarks justify the statement made at the beginning of this 
section, namely that the spin network simulates intrinsically time evolution 
(without resorting to ad hoc external Hamiltonians as happens for instance 
in the Universal Quantum Simulator model proposed in ^0]). Moreover, 
by changing the program V (and the algorithm A) the machine is able to 
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select different types of dynamical behaviors, and thus to simulate complex 
poly-local many angular momenta interactions modeled as binary couplings 
and Wigner rotations (see the assumptions discussed in Section 2). What we 
mean in particular is that different types of evolutions in 0„ can be grouped 
into "computing classes" based on the choice of gates that each program 
has to employ. Then a program V (defined in (j59|) and associated with a 
collection of directed paths as in (jUUI) or, equivalently, with decompositions 
into a number of gates as in ()63p) may be based on either M-gates alone, or 
j-gates alone, or some fixed sequence of M and j-gates. 

♦ An M- computing class contains programs which employ only M-gates 
at each step in their associated directed paths. The binary bracketing struc- 
ture of the computational Hilbert spaces described in Section 2 is not in- 
volved, and it is not difficult to realize that such kind of computation, when 
applied to N ^-spins, reproduces the usual Boolean quantum circuit (c/r. 
the end of Section 3.2 and Appendices Bl and B2).A 

♦ A j-computing class includes programs which employ only j-gates at 
each computational step. This class is particularly interesting since it shares 
many features with suitable types of discretized field theories (the so-called 
state sum models, to be discussed in Section 6) as we already noticed in 
[TTj . Now the combinatorial structure of Rotation graphs becomes prominent 
owing to the existence of an one-to-one correspondence between allowed 
elementary operations and the edge set E of 0„, for each n [cfr. pUj) and 
more generally Appendix A2). 

In the present framework it is convenient to switch back to notations used 
in the first part of Section 4.1. Then states will be labeled again by |b)„ as 
in (j43p (dropping JM ) and a program V is represented formally as 

■p 

Wsnj [bin — ^ bout] : |bin)n ^ |bout)n, (70) 

where, as before, |bin)n is fixed and |bout)n is an accepted state. The set 
of operators Hp in has been replaced by the unitary operators U^nj 
introduced in (fTUI) and already used in (j^ . The collection of directed paths 
associated with (fTUj) is defined as in (p?|) and a particular path of length L 
corresponding to a time-ordered sequence of states is represented as 

I bin )n = \bo)n ^ \bl)n ^ ^ \ bs)n ^ ^ I bL )n = I bout )n (71) 
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Each arrow corresponds now to the first type of operation in (j62p . namely a 
Racah transform (possibly up to weights/phases). When one particular path 
is chosen we would recover expressions similar to referring to a 

computation process based on a program V performing an algorithm A. 

However, in the j-computing class one may address other types of prob- 
lems, namely: selected two states in say |bin)„ and |bout)n5 consider all 
possible V{A) that compute |bout)n as the result of the application of some 

W3„,[bin ^ bout] (72) 

to |bin)n- The functional on 0„ which takes care of such multiple choices is 
a "path sum" (a discretized Feynman's path integral) which may be written 

as 

Z [bin, bout] = 2^ (bout I W3„j[bin ' bout] | bin )„ , (73) 

V{A) 

where the summation is over all paths with fixed endpoints and W-pi^^) is a 
weight to be assigned to each path. 

Notice that if we should give the same weight, say PVp(^) = 1 to each path, 
then the results on equi-probability amplitudes collected in the theorem of 
Section 4.1 ensure us that the functional ()73p is a combinatorial invariant, 
namely it is actually independent of the particular path connecting bin and 
bout- On the other hand, if we insert non trivial weights in ()73|). we may natu- 
rally address questions about most efficient algorithms and time complexity. 
For instance we could weigh paths with the inverse of their lengths L{V{A))] 
then the minimum-length path (the optimal algorithm) will be dynamically 
singled out in the path sum. As a matter of fact, even such simple example 
turns out to be highly non trivial owing to the combinatorial complexity of 
0„. We shall come back on such issues in the following Section 4.3. A 

♦ An altered j-computing class is a modification of the j-class obtained 
on applying just one Wigner rotation (M-gate) to the input state. A typ- 
ical directed path in this class can be represented by mixing our previous 
notations (see (fHH|) . (fHH|l . and (ITTll'l to sret 

|Oin)n = |bo;6',0; J)n^ |bo;6'',0'; J)„ ^ |bi )„ ^ |b2 )« ^ ■ • ■ ^ 

^ |b.)n ^ > \bL+l]e\<P']J)n = |Dout)n (74) 
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where the first arrow represents a Wigner rotation 2)"^(a/?7) and all the others 
are alterations of the binary bracketing structure while keeping the angular 
dependence of the states fixed. Such computing class seems sufficiently gen- 
eral and shares some features with the j-class for what concerns the path 
sum interpretation. The counterpart of the path sum functional (|7Hjl reads 

Wp(_A){bL+l\K3nj[bo ^ bL+l]\bo)n 

V{A) 

{bo;e',(j)';J\V'{aP^)\bo;e,(j);J)^, (75) 

from which we see in particular that combinatorial invariance is broken even 
if we assign to each path the same weight. A 

♦ An alternating computing class includes programs which alternates M 
and j-gates and the length of each of the associated directed paths is an 
even number L = 2i. This class is quite general with respect to capability 
of simulating real physical systems and includes all the former computing 
categories since anyone of the elementary gates could eventually be realized 
by an identity transformation. A 



4.3 Computational complexity 

In view of the role that binary coupling trees play in our model for quantum 
computation we define an encoding map 

K^(b) ^ Tib) (76) 

where, as in (15), the shorthand notation b stands for the string of quan- 
tum numbers ([ji, j2, • • • jn+i]''; ^i, • • • 5 ^n-i); T(b) is the coupling 
tree uniquely associated with the computational Hilbert space for given 
J and n (see Appendix Al). This coding is intrinsically quantum, namely 
T(b) in (96) is not a simple device in which classical information can be stored 
(as happens for instance when search trees are considered). The quantum 
behaviour of T(b) emerges in particular when we look at the nature of the 
internal labelings. Recall from Section 2 that an intermediate angular mo- 
mentum, say Ki = Ji + J2, has eigenvalue ki ranging between \ji — J2I and 
ji + J2 and thus such quantum trees - even when equipped with definite values 
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of the incoming quantum numbers chosen in {0, |, 1, |, . . .} - take care con- 
sistently of the range of different values that may be assigned to the internal 
nodes. 

It is worth noting that the binary bracketing notations introduced in (12) 
and (14) of Section 2 can be interpreted as the quantum counterpart of the 
"word construction" outlined in Remark A.l of Appendix Al for a generic 
binary operation. To formalise this observation we introduce a new map - 
equivalent to ()76|) from a quantum mechanical point of view - which encodes 
information carried by the Hilbert space Ti^ into a quantum word, namely a 
string of quantum numbers plus parenthesization. Denote formally this map 
by 

n^w ^ ([ji,j2,...jn+i]'; fc^fc2^...,fctl), (77) 

since a more explicit form as in (14) would force us to write down just 
one particular type of binary coupling. In our opinion a closer inspection 
of the encoding map (f77j) (instead of ()76|) which will be exploited in the 
following) could represent a promising starting point to establish a truly 
quantum Formal Theory including languages, grammars, Godel numberings 
and related automaton models [TT j. 

Coming back to the encoding map ()76|). and referring to the topological 
transformations on binary coupling trees discussed in Appendix A2, we easily 
recognise that a Racah transform TZ defined in (21) is encoded in a rotation 
(more precisely, the explicit expression given in (23) is encoded into the 
operation depicted at the bottom of Fig. 19). On the other hand, a phase 
transform $ defined in (22) turns out to be encoded into a twist (compare 
e.g. the explicit expression (24) with the twist depicted at the bottom of Fig. 
20). Consequently, 

Racah transform TZ — > Rotation 
Phase transform $ — > Twist (78) 

represent encoding maps associating the two types of unitary elementary j- 
gates introduced in Section 3.1 with basic topological moves on quantum 
trees. The role of these two sets of operations is specular also for what 
concerns composition, since the Biedenharn-Louck Theorem (Section 3.1) 
has its counterpart in the fact that any pair of binary coupling trees can be 
connected by a sequence of rotations and twists (see Appendix A2). 
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The quantum encoding maps (j76|) . (|78p make manifest that combina- 
torics of (Twist)-Rotation graphs (Appendix A2) and of the computational 
space of the simulator (Section 4.1) share identical features, at least as far 
as j-computing classes (Section 4.2) are implemented. This crucial remark 
justifies the fact that we may speak about combinatorial and computational 
complexity questions by employing a common language and concepts. Note 
however that these similarities can be exploited only to some extent since 
combinatorial complexity of graphs is usually addressed in a classical in- 
formation theory context (see Appendix A3). In the remaining part of this 
section we shall illustrate in brief the computational potentialities of the spin 
network simulator. Going beyond questions in number theory, we argue that 
our model is suitable to deal with #P ("hard enumerative/combinatorial") 
problems more efficiently than any classical machine. 

For what concerns space complexity capacity of the spin network sim- 
ulator in the sense of capability of storing information, we realize that it 
behaves as the cardinality of the (Twist) -Rot at ion graphs (c/r. ()114|) . ()115|) . 

Table 2 and Table 3 of Appendix Al). When the number n of incom- 
ing angular momenta grows, the number of states which becomes accessible 
for computation increase at least exponentially. To quantify these asymp- 
totic growth, consider first the Catalan numbers (96) which represent some 
sort of lower bound for the various enumerations of rooted labeled binary 
trees shown in Table 3 of Appendix Al. They have the following asymptotic 
expansion for n ^ oo 

gnln4 [11 1 

^" ^ ^^n{n + if r"8;^+128;^ + --7 ^^^^ 

where In = logg. On the other hand, we may estimate approximatively rates 
of grows of the double factorial Dn in ()116|) and of the quadruple factorial 
Cn in ()115|1 by using Stirling formula. We get 

D„ ?s n " exp {n In 2 - n} (80) 

C'„ n" exp{2nln2 - n} (81) 

where the subleading terms are decreasing (increasing) exponentials, respec- 
tively. Thus the case considered in Section 4.1, namely the computational 
space modeled on the Rotation graph C5„, turns out to exhibit for large n 
space complexity of factorial class ~ n! as in (jHO)). 
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To address analogously time complexity, we need first the notion of "input 
length " (for an instance of some given problem) which turns out to be related 
with the encoding scheme employed. By using the map ()76|) we may say 
that a typical input length is the number of symbols required to specify a 
(quantum) labeled tree, namely (2n + l) (terminal nodes, intermediate nodes 
and the root). Thus it seems natural to assume such a number as a typical 
measure of the size of the input (note however that in order to specify one 
particular quantum state we should choice also a value for the total magnetic 
number M). Although the input length is linear in the number of symbols, 
the quantum nature of the computational space is reflected by the fact that 
the size of the configuration space accessible for computation grows factorially 
with n as discussed in the previous remark. 

With these premises and by exploiting the estimate on the Diameter of 
the rotation graph &n given in ()128p of Appendix A3, we assert that the 
(time) complexity function for any possible algorithm running over the spin 
network simulator can be expected to be polynomially bounded as a function 
of the input size n. An effective discussion of both space and time complexity, 
however, requires of course reference to a specific algorithm, which in turn 
can be formulated only provided the necessary encoding scheme is defined. 
Work is in progress along these lines [T7j . 



5 Semiclassical simulator 

and SU(2) state sum models 

According to the Bohr correspondence principle, classical concepts become 
increasingly valid in the regime where quantum numbers are large. In han- 
dling with angular momenta variables measured in units of h, the classical 
limit h ^ implies that, for finite angular momenta, both the j-quantum 
numbers and the magnetic ones are much bigger than one. For what concerns 
pure angular momentum states - and in particular the computational Hilbert 
spaces introduced in Section 2 and involved in dynamical processing as illus- 
trated in Section 4.2 - when approaching classical limit all the components 
of the operators {Ji{i = 1, 2, . . . , n -f- 1), J} are confined to narrower ranges 
around specific values. Thus geometrical concepts typical of the semiclassical 
vector model arise naturally and the corresponding physical quantities have 
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to be thought as averaged out. As we shall see below, angular momentum 
functions such as Racah transforms and Wigner rotation matrices admit well 
defined asymptotic limits, whose absolute squares (probabilities) correspond 
to classical limits of the related physical quantities. 

With these preliminary remarks, and on the basis of JHl and Topic 
9 (in which a self contained discussion of the various asymptotics is given, 
together with the list of original references), we are going to set up a semi- 
classical counterpart of the spin network simulator which represents, to our 
knowledge, the first explicit example of a quantum circuit model mapped 
onto a (classical) probabilistic automaton scheme. 

Let us focus for the moment on probabilities, namely on asymptotic ex- 
pansions of absolute squares of transition amplitudes. Consider an elemen- 
tary j-gate, namely a Racah transform expressed in terms of a 6j symbol as 
in (j2ni)- When all the six angular momenta in the 6j become 3> 1 in units, 
the square of the symbol has the limiting value given by the Wigner formula 

a b 



where V is the Euclidean volume of the tetrahedron formed by the six angular 
momentum vectors whose lengths are the arguments of the coefficient (V^^ 
can be computed from (a, 6, c, e, /) by using the Cayley determinant). This 
result can be exploited to find the probability of measuring a coupling scheme 
(a(6c)e)/ having prepared the system in the scheme {{ab)dc)f . Denoting by 
P(e) this probability and using and (jH^ we get 



^ (2e + l)(2d+l)^ 

^ ' 12nV ^ ^ 

Coming to elementary M-gates, and in particular to the expression 
for a W-rotation matrix in terms of Euler angles, we see that 

I^m'mMt)!' = MM'M(/5)r (84) 

and this quantity is symmetric in M' and M. Following the step illustrated 
in [in] (Topic 9, Section 2) we limit ourselves to analyse the case of an 
input quantum state characterised in the original reference frame by a total 
angular momentum J maximally oriented along the z-axis, namely M = J. 
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the probability that the angular momentum projection along the rotated 
2;'-axis has the value M' = M is given by 



P{M) = [dl,^{f3)r = ( , ,J cos^ sin^^ . (85) 



J - MJ \ 2 J \ 2_ 

In the classical limit the most probable value of M is distributed around the 
classical value Mq = J cos (3 with some fixed value of the classical probability 
P{Mo). Then it can be shown that the limiting value of (jH3jl for J ^ 1 and 

J ± Mo > (M - Mo) reads 



P(M) ~ P(Mo) exp <^ - - . ^ ° } . (86) 




Thus the probability for M is a Gaussian distribution around the classical 
value Mo and the dispersion in the variable (M — Mo), a/J sin/5, is of the 
order \/j by the assumption made above. The simplicity of this result is due 
to the choice of | J J > as original states. More general types of classical limits 
and asymptotics are discussed in [13j (Topic 9, Section 10) and collected also 
in HH. 

Without entering into further technical details, the rationale underlying 
our approach should have become clear: anyone of the computing classes of 
the quantum simulator introduced in Section 4.2 as finite sequences of el- 
ementary unitary j- and/or M-gates can be mapped onto a corresponding 
(classical) non-deterministic circuit-type computing process based on com- 
binations of (jH^ and/or (jH^ . We argue that such a semiclassical model 
could be able to simulate physical systems made up by many interacting 
constituents such as polyatomic molecules described by pure states of some 
suitable angular momentum-type variables for high values of q-numbers. 

In order to complete our semiclassical picture we have to include the treat- 
ment of asymptotic limits of the transition amplitudes associated with the 
elementary gates employed in the quantum model. As we shall see, the result- 
ing setting is closely related to SU{2) state sum models introduced in a com- 
pletely different context, namely topological quantum field theories (TQFT) 
and Euclidean quantum gravity defined on triangulated 3-dimensional man- 
ifolds (see 1221, foi' extended reviews on such topics). 
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The key point is the interpretation of the Ponzano-Regge asymptotic 
formula for the 6j symbol which reads ^H] 



where the limit is taken for all entries ^ 1 (recall that ^ = 1) and ir = jr^^/'^ 
with {jr} = {a,b,c,d,e, f}. V is the Euclidean volume of the tetrahedron 
with edges of lengths {ir} (note the shift j ^ j + 1/2 with respect to the 
variables employed in calculating the volume in (jHH)) and finally Or is the 
angle between the outer normals to the faces which share the edge ir- 

• From a purely quantum mechanical point of view, the probability am- 
plitude (jHZj) has the form of a semiclassical (wave) function since the 
factor l/v247rK is slowly varying with respect to the spin variables 
while the exponential is a rapidly oscillating dynamical phase. Such 
behavior complies with the fact that the square of the modulus of 
the asymptotic (jH7|) reproduces Wigner's expression (jH^ . Moreover, 
according to Feynman path sum interpretation of quantum mechan- 
ics, the argument of the exponential represents a classical action, and 
indeed it can be read as ^ p g for pairs {p, q) of canonical variables 
(angular momenta and conjugate angles). 

• There exists another intriguing physical interpretation of (|87|) if we 
recognise that the expression in the exponential represents the classical 
Regge action [20] ^ namely the discretized version of Einstein-Hilbert 
action of General Relativity - for the tetrahedron associated with the 
6j symbol in the asymptotic regime. 

In Regge's approach the edge lengths of a triangulated spacetime are 
taken as discrete counterparts of the metric tensor appearing in the 
usual action for gravity and angular variables (deficit angles) are related 
to the scalar curvature obtained from the Riemann tensor. Strictly 
speaking, a "triangulated spacetime" is a piecewise linear (PL) mani- 
fold of dimension D dissected into simplices, namely triangles va. D = 
2, tetrahedra in D = 3, 4-simplices in D = 4 and so on. Inside 
each simplex either an Euclidean or a Minkowskian metric can be as- 
signed: accordingly, spacetime manifolds obtained by gluing together 
D-dimensional simplices acquire an overall PL metric of Riemannian 
or Lorentzian signature. 



39 



The Regge Calculus formalism became in the early 80's the start- 
ing point for a novel approach to quantization of General Relativity 
known as Simplicial Quantum Gravity (see the review [31] and refer- 
ences therein). The quantization procedure most commonly adopted 
is the Euclidean path sum approach, namely the discretized version of 
Hartle-Hawking path integral describing D-dimensional, locally Eu- 
clidean geometries undergoing "quantum fluctuations", possibly with 
the constraint of keeping some {D — 1) -dimensional boundaries fixed. 

Coming back to the interpretation of (j87|) . we conclude that it repre- 
sents the semiclassical functional - to be intended as the semiclassical 
limit of a sum over all quantum fluctuations - associated with a very 
simple 3-dimensional "spacetime", the Euclidean tetrahedron. 

On the basis of the remark above, we pass to describe in brief the Ponzano- 
Regge state sum model representing the (quantized) partition function of 
simplicial Euclidean 3-gravity. Denote by 

T\j) (88) 

a particular triangulation of a closed 3-dimensional PL manifold M.^ (of fixed 
topology) obtained by assigning SU{2) "spin variables" {]} to the edges of 
. The assignment must satisfy a number of conditions which can be more 
easily illustrated if we introduce the state functional associated with T^{j), 
namely 



z[r3(j)--M^^]=A(^)~^°n(-i)''"^^n<^M? ? ? ^89) 

where A'^o, A^i, denote the number of vertices, edges and tetrahedra in 
r3(j), A(L) = 4LV3C (C an arbitrary constant), = {2jA + 1) are the 
dimensions of irreducible representations of SU{2) which weigh the edges, 

(pB = (—l)'^p=^^p and 1^ ^ ^1 are SU{2) 6j symbols to be associated 

with the tetrahedra of the triangulation. The Ponzano-Regge state sum is 
obtained by summing over triangulations corresponding to all assignments 
of spin variables {j} bounded by the cut-off L, namely 

Zpn[M'] = hm V Z[T3(j)^A<3.^]^ (90) 

{j}<L 
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where we formally remove the cut-off by taking the limit in front of the sum. 
As already noted in [I^, the above state sum is a topological invariant owing 
to the fact that its value is actually independent of the particular triangu- 
lation, namely does not change under suitable topological transformations 
(the bistellar moves). These moves are expressed algebraically in terms of 
the Biedenharn-Elliott identity (j46p -representing the moves (2 tetrahedra) 
^ (3 tetrahedra)- and of both the Biedenharn-Elliott identity and the or- 
thogonality conditions ()48|) for 6j symbols, which represent the barycentric 
move together its inverse, namely (1 tetrahedra) ^ (4 tetrahedra). 

The state sum (jUUj) (and, more generally, geometric partition functions of 
this type built up in any dimension D [3.5 J resembles the functional (|73p in- 
troduced in dealing with the j-computing class of the spin network simulator 
(Section 4.2). This is due to the fact that the amplitude of the 3nj symbol 
in (fTHj) can be factorized according to the general prescription (jUHjl into sums 
over intermediate angular momenta of products of 6j symbols weighted by 
suitable factors and phases (c/r. fOHj) . ((211) )■ These two partition functions 
share the property of being combinatorially invariant under topological moves 
expressed in terms of algebraic identities of the 6j symbols. However, on the 
one hand, the Racah identity does not appear in the Ponzano-Regge 
framework since it would correspond to a topological transformation (1 tetra- 
hedron) ^ (2 tetrahedra) which is forbidden in the PL category. On the other 
hand, in the spin network framework it is not required a priori that the 6j 
symbols match together to give rise to a triangulation of a 3-dimensional 
manifold. Moreover, if we fix n to get a specific computational space mod- 
eled on the graph (S„, we would not catch in (J5„ all possible triangulations of 
a given PL 3-manifold. Although we may be tempted to claim that the spin 
network is able to simulate 5'f/(2)-coloured 3-dimensional quantum gravity, 
we should bear in mind that we are actually dealing with a graphical device 
which encodes all types of 3nj symbols for any fixed n [22j. In this per- 
spective it is interesting to recall that Ponzano and Regge themselves 
noted that the topology of a 9j symbol corresponds to the real projective 
space MP^, in the same sense that the 6j has the topology of the 2-sphere 
bounding the tetrahedron. Indeed any particular type of 3nj symbol may 
be associated with a closed, not necessarily oriented, surface representing 
the boundary of a 3-dimensional polyhedron |36j obtained by duality from 
the graphical representations introduced in ^21- In this sense the claim that 
the simulator can simulate some classes of extended triangulated objects is 
certainly true in dimension 2. 
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Remark 5.1. In |SZj a regularized version of -based on representa- 
tion theory of a quantum deformation of the group SU (2)- was proposed and 
shown to be a well-defined (finite) quantum invariant for closed 3-manifolds. 
Its expression reads 

ZTv[M^;q] = Yi ^""^^ n^^ n 

{j} A=l B=l 

where the summation is over all colourings {j} labeling highest weight ir- 
reps of SU{2)q {q = exp{27rVr}, with {j = 0, 1/2, l...,r - 1}), = 
(— l)^-^-4[2j^ _l_ x]^ where []q denote a quantum integer, w = 2r/(g — g^^)^ 
and I \b represents the q-6j symbol whose entries are the angular momenta 
je , i = 1, . . . , 6 associated with tetrahedron B. In [38j the invariant ()9ip is 
shown to equal the square of the modulus of the Reshetikhin-Turaev invari- 
ant, which in turn represents the Chern-Simons partition function written 
for a closed oriented manifold equipped with a surgery presentation. We 
may write schematically 

ZTv[M';q] ^ \ Zcs[M^;k]\\ (92) 

where the level k = 2(r — 1) of the Chern-Simons functional is related to 
the deformation parameter q and also to the cosmological constant of the 
underlying Euclidean gravity model (c/r. |22 (Ch 7) and ^\ for reviews on 
Chern-Simons theory and its relations with 3-dimensional gravity) . 

As we shall see in the following section, functors derived from SU{2) 
Chern-Simons theory are the basic ingredients for implementing computa- 
tion in the topological approach It will be shown that the spin network 
dynamics based on j-gates can be mapped into the functorial approach in a 
way that resembles the correspondence (|92j). ▲ 



Jl J2 J3 

ii ib je 



(91) 



6 Spin network and topological 
quantum computation 

We begin this section by introducing some basic ingredients of Chern-Simons- 
type Topological Quantum Field Theories (CS TQFTs) in order to deal with 
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the topological approach to quantum computation. Our presentation will 
be necessarily sketchy, and we refer the reader to jlHI; IH]; IHI for general 
reviews on TQFTs, while the 3-dimensional CS case is extensively addressed 

in [sg, El, M- 

TQFTs are particular types of gauge theories, namely field theories quan- 
tized through the (Euclidean) path integral presciption starting from a clas- 
sical Yang-Mills action defined on a suitable D-dimensional space(time). 
TQFTs are characterized by observables (correlation functions) which de- 
pend only on the global features of the space on which these theories live, 
namely they are independent of any metric which may be used to define 
the underlying classical theory. The geometrical and topological generating 
functional and correlation functions of such theories are computable by stan- 
dard techniques in quantum field theory and provide novel representations of 
certain global invariants (for D-manifolds and/or for particular submanifols 
embedded in the ambient space) which are of prime interest in mathemat- 
ics. In the 3-dimensional case, theories based on Chern-Simons-type actions 
(see below) have been shown to incorporate significant generalizations of pre- 
viously known invariants for both 3-manifolds (Witten-Reshetikhin-Turaev 
invariant) and knots/links (Jones polynomial). In particular, the Jones poly- 
nomials can be obtained as correlation functions of Wilson line operators 
along closed loops in the CS framework jUj. While these mathematical ad- 
vances are self-evident, CS theory also provides a unifying 3-dimensional 
viewpoint for 2-dimensional Conformal Field Theory as well as new results 
on 3-dimensional quantum gravity. 

Since TQFTs are quite generally soluble, they could provide a testing ground 
for new approaches to the quantum theory of fields. It has been conjectured 
that TQFTs may represent different "phases" -in which general covariance 
is unbroken- of their more conventional counterparts. 

Denote by Si and S2 a pair of 2-dimensional manifolds and by Ai^ a 
generic 3-dimensional manifold with boundary dAi^ = S1US2 (all manifolds 
here are compact, smooth and oriented). A unitary 3-dimensional quantum 
field theory corresponds to the assignment of 

i) finite dimensional Hilbert spaces (endowed with non-degenerate bilinear 
forms) and Ht;2 to Si and S2, respectively; 
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ii) a map ("functor ") connecting such Hilbert spaces 

He. ^ H., (93) 

where Ai^ is a manifold which interpolates between Si (incoming boundary) 
and S2 (outgoing boundary). Without entering into details concerning a few 
more axioms (diffeomorphism invariance, factorisation etc.) we just recall 
that unitarity implies that 

iii) if S denotes the surface S with the opposite orientation, then Hf; = Ti"^, 
where * stands for complex conjugation; 

iv) the functors (jHSl) are unitary and Z[7W^] = Z*[7W^], where denote 
the manifold with the opposite orientation. 

In Chern-Simons theory the functor (j93|) is the partition function Zqs 
associated with the classical action 

Scs{A) = [ (AdA+-AAAAA] (94) 

written for simplicity for a closed manifold Ai^, dAi'^ = 0. Here A is an 
SU (2)-connection, namely a 1-form on the principal S't/(2)-bundle over A4^, 
d is the exterior differential and A is the wedge product of differential forms. 
The partition function is obtained functionally by integrating the exponential 
of the classical action ()94j) over the space of all SU{2) connections. We 
formally write 

Zcs[M';k] = I [DA] exp | Scs (A) | (95) 

where the coupling constant k (the level of the theory) must be an inte- 
ger. It can be shown that the partition function ()95|1 actually represents a 
topological invariant for closed 3-manifolds, related in turn to the Turaev- 
Viro invariant by the correspondence (where |Zcs'[A1^]p stands for 

The extension of (j95j) to the case dAi^ 7^ requires modifications of the 
classical action by suitable (Wess-Zumino-type) boundary terms [HU] . 
In view of applications in a computational context it is sufficient to note that 
for a boundary component S, Z(75[A^^](E) can be realized as restriction 

to invariant subspace W C Tis of transformations of the form JJ^fi'j on 

i 

the subspace of computational states of a quantum computer, where the Qi, 
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interpreted as gates in a quantum circuit scheme, can be written as "words " 
in the standard generators (Dehn's twists) of the Mapping Class Group of 
the surface S. Also in the computational framework the "observables" of the 
theory turn out to be represented by Wilson loops as discussed in [45j and in 
the following paragraph. The explicit expression of Wilson loops operators 
yVk{K), namely holonomies of the connection 1-form evaluated on closed 
curves K in read 



where hoi K = V exp / A, V is the path ordering, and A is now thought of 

Jk 

as the connection over the ^-bundle of Lie algebra-valued 1-forms tangent 
to M.^ {Q being the gauge group). 

It is worth to recall that the evaluation of Jones polynomials - appearing in 
expressions like - was shown to be computationally #P jlHj (namely 
essentially the enumerative equivalent of NP-complete [TF]). 



6.1 Holonomic Quantum Computation 

Holonomic Quantum Computation (HQC) is an all-geometrical approach 
to quantum information processing. In the HQC strategy information is 
encoded in degenerate eigenspaces of a parametric family of Hamiltonians. 
The computational network of unitary quantum gates is realized by driving 
adiabatically the Hamiltonian parameters along loops in a control manifold. 
By properly designing such loops the non-trivial curvature of the underlying 
bundle geometry gives rise to unitary transformations, i.e. holonomies that 
implement the desired unitary transformations. Conditions necessary for uni- 
versal QC are stated in terms of the curvature associated to the non-Abelian 
gauge potential over the control manifold. In view of their geometrical na- 
ture the holonomic gates are robust against several kind of perturbations 
and imperfections. This, along with the adiabatic fashion in which gates 
are operated, makes in principle HQC an appealing way towards universal 
fault-tolerant QC. 

HQC as introduced in 0, [IHl, is based on a novel gauge-theoretic frame- 
work in which one is supposed to be able to control a set of parameters 



Wk{K)= f [DA]e^^^^(^)Tr(holir) / / [DA] 



(96) 
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A G £, on which depends an iso-degenerate family T of quantum Hamilto- 
nians \H{X)\. Information is encoded in a //-dimensional eigenspace C of a 
specific H{Xo) G JF. Universal QC |48| over C can be then obtained by adia- 
batically driving the control parameters along suitable loops 7 rooted at Aq. 
The key physical ingredient is provided by the appearance in such quantum 
evolutions of non-Abelian geometrical contributions [IH] ^ > 1) 

given by holonomies associated with a gauge potential A valued in the al- 
gebra of f/(z/)[5TJ", In other words quantum computation in the HQC 
approach is nothing but the parallel transport of states in C realized by the 
connection A. Therefore the computational power in the HQC approach 
relies on the non-triviality of the geometry of the bundle of eigenspaces of 
JF over the manifold of control parameters, C: in this sense HQC is fully 
geometrical. It is worth observing that the computational subspace C can 
be thought of as the lowest-energy manifold of a highly symmetric quantum 
system; from this point of view HQC is a kind of ground-state computation. 
This last remark points out the potential existence of a fault-tolerant |^ 
feature of HQC due to energy gaps and even spontaneous relaxation mecha- 
nisms. Further fault-tolerant characteristics of HQC are related to the fact 
that the holonomies U-f realizing quantum computations typically turn out 
to depend just on the areas of the surfaces that the generating loops 7 span 
on certain 2-dimensional submanifolds. When this area is given one can 
consider even very large, i.e. "far" from the identity deformations of 7, but 
as long as they are area-preserving no errors are induced. Moreover as far 
as the adiabaticity condition holds, does not depend on the rate at which 
the control loops are driven. Hence, even with respect the issue of timing, 
HQC can be expected to be robust. 

The evolution of the quantum system is thought of as actively driven by 
the parameters A, over which the experimenter is assumed to have direct 
access and control, being able to drive by a dynamical control process the 
parameter configuration A G £ through a path 7 : [0, T] — >• C Hence, a 
one-parameter (time-dependent) family 

= {H{t) = o 7(t)]|t G [0,T]} C ^ , (97) 

is defined for all $ : £ 1-^ U{Af) , Af = dimC, $ being a smooth mapping, 
and dim ([/(A/")) = N"^. The quantum evolution associated to the family 
J-'^ is described by the time-dependent Schrodinger equation idt\ip{t)) = 
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H{t)\'tl){t)) and hence it has the operator form 

C/7 = T exp I ^ dtH{t) I e U{Af) , (98) 

where T denotes chronological ordering. The above time-dependent quan- 
tum evolution, for a given map depends in general on the path 7 and not 
just on the curve 7([0, T]), namely the image of 7 in the control manifold. In 
other words the unitary transformation f/^ contains a dynamical as well as a 
geometrical contribution, the former depends on the rate at which 7([0,T]) 
is traveled along whereas the latter depends merely on the geometrical char- 
acteristics of the curve. 

Non-Abelian holonomies are a natural generalization of the Abelian Berry 
phases. The basic assumption is that JF is an iso-degenerate Hamiltonian 
family, i.e. all the elements of JF have the same degeneracy structure. When 
the system control parameters are driven adiabatically, slowly with respect 
to any time-scale associated to the system dynamics, along a loop 7 in £ any 
initially prepared state \ipin) £ 'H will be mapped after the period T onto the 
state 

R 

l^out) = C/7 \M , = e''^^r^,(7) , (99) 

e=i 

where (f)^— I dTei{\T-) is the dynamical phase {£e{X) denoting the degen- 

erate Hamiltonian eigenvalues) whereas the matrices r^^(7) represent the 
geometrical contributions. They are unitary mappings of Hg, onto itself and 
they can be expressed by the following path ordered integrals 

TA,{n) ^V(fAe e U{ne) , £ = 1, . . . , , (100) 

denoting the number of degenerate states of energy ei, and R the number 
of degeneracies. These are the holonomies associated with the loop 7, and 
the adiabatic connection forms A^. The latter have an explicit matrix form 
given by 

Ae = E^^./^dA, , {AeX" = m>^)\^\i^i'W) , (101) 
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with {^fj.} the local coordinates on C, and Hi = span{?/'"(A))}"^^, 
ip"{X)) denoting the eigenstates of H{\) corresponding to eigenvalue ^^(A). 
The connection forms Ai are nothing but the non-Abelian gauge potentials 
enabling the parallel transport over C of vectors of the fiber Tii. 

The combinatorial setting of the spin network simulator, and in particular 
of its fiber space structure over (see Remark 4.1 at the end of Section 4.1 
and 29 ), appears to provide the natural structure for a discrete-time im- 
plementation of HQC, at least when the Hamiltonian defined in (I69|) exhibit 
the required degeneracy. 

6.2 Combinatorial setting of Topological 
Quantum Computation 

The approach to quantum computation of Freedman and collaborators j3] 
is based on an extension of Chern-Simons functor at level = 3 (CSS) 
to particular types of 2-dimensional boundary objects, namely closed disks 
with some marked points. Denote by {D^, 3 pts) a closed disk -namely a set 
topologically equivalent to the standard 2-disk {(x, G + < 1}- 

with three points lying in the interior of D^. These points, together with 
the boundary dD^ of the disk, are marked by four labels {a, b, c, d} chosen in 
the set {0, 1, 2, . . . k}, which reduces to {0, 1, 2, 3} in the present case k = 3. 
Notice that the convention adopted by the authors of ^ is different from 
ours: the level k = 3 here corresponds to their k = 5. 

The disk considered so far and depicted in Fig. 3 (left) is the support of the 
so-called "topological qubit ". More precisely, a topological qubit corresponds 
to the image of the usual 1-qubit space into the Hilbert space TCcs3 (-D^, 3 
pts) derived from CS theory. The inclusion 

withi(C2) = 1/(D^3pts) (102) 

induces on the label set {a, b, c, d} the identifications a = b = c = 1 and 
d = 0, namely a mapping onto the binary digits {0, 1} (see Fig. 3, right). 

Figure 3 
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In the present framework 1-qubit gates are implemented by considering - 
instead of the full Chern-Simons functor (|93|) - a unitary action of the braid 
group i?3 on the strands generated in the ambient space by "evolution" of 
the marked points. 

Figure 4 

If 

^3 X (L)^3pts) ^ (L'^3pts) (103) 

denotes the action of the braid group -B3 on the topological support (see Fig. 
4) then we formally write the induced unitary functor between Hilbert spaces 
as 

Vi(D2^3pts) V2(^^3pts). (104) 

A configuration of iV qubits is supported by a 2-disk with 3iV marked points 
(D^, 3A^ pts) and the associated Hilbert space is V^D"^, 3N pts), namely the 
image of (C^)®-^ in TCcs5 {D^, 3A^ pts), where the image is defined as in ()102|) . 
An A^-qubit gate is represented by a unitary map 

Vi{D\3Npts) V2(^^3A^pts), (105) 

where B^j^ is the braid group acting on the 3A^ strands generated by the 
marked points. 

We may summarise the main results of the series of papers [4 as fol- 
lows. It is shown that their CS'5-functor (and also any other C5'/c-functor, 
k > 3, k ^ 4) is universal for quantum computation and in particular that 1 
and 2-qubits topological gates are sufficient to reconstruct all other gates of 
the type ()105|) . The resulting model is polynomially equivalent to the usual 
Boolean quantum circuit and, conversely, it is actually shown that there 
exists (at least) a class of TQFTs which can be simulated on a quantum 
machine (c/r. the diagram at the end of Section 2). There are however some 
open problems in this approach, mainly due to the fact that it is difficult to 
"localise" topological objects (such as the marked points on disks) in order 
to provide local Hamiltonian operators. 
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Remark 6.1. It may be useful to recall some basic properties of the 
Artin braid group -B„ (see for instance |S2|)- Bn has n generators, denoted 
for the moment by {ai, (72, ... , a„}, which satisfy the relations 

(Ti(Tj =(Tjai if \i- j\ > 1 
(Ti dj+i (Tj = (Ti=i (Ji (Ji+i (i = 1,2, . . . ,n). (106) 

This group acts naturally on topological sets of n disjoint strands - running 
downward and labeled from left to right - in the sense that each generator ai 
corresponds to a crossing of two contiguous strands labeled by i and (z + 1), 
respectively (if stands for the crossing of the i-th strand over the (i + 1)- 
th one, then represents the inverse operation and CTj = (J~^ai = 
Identity). By slightly changing notations, denote by Rij the (over) crossing 
operation acting on two strands the endpoints of which are labeled by i and 
j. Then the second relation in (|106|) can be recasted into the form 

-R12-R13-R23 = -R23 -R13 -R12 (107) 

and represented pictorially as in Fig. 5, where operations are ordered down- 
ward. Note that this picture can be viewed as a portion of an n-strands con- 
figuration (and thus {1,2,3} may actually represent labels attached to any 
triad of contiguous strands) since the first relation in ()106|) ensures that other 
kinds of crossing are trivial. The relation ()107|) (or, alternatively, ()106p ) is the 
algebraic Yang-Baxter equation which characterises the algebraic structure 
of a number of models in statistical mechanics and field theory. 

Figure 5 

The Artin braid group -which arises naturally in the CS-framework and 
consequently in the topological setting for quantum computation described 
above- plays a crucial role also in the approach recently proposed in 
where unitary representations of B4 are shown to be universal gates for 
Boolean quantum computation. We shall come back on this point at the 
end of this section, after the analysis of the algebraic structure underlying 
the spin network model. A 

We come now to describe the algebraic content of the theory underlying 
the spin network simulator by showing explicitly how binary coupling trees 
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can be embedded into "combinatorial " 2-disks with marked points whose as- 
sociated transformations satisfy the Racah and the Biedenharn-Elhott iden- 
tities (jmi) and (jlTj) (instead of the Yang-Baxter equation p07p characterising 
the standard topological approach). 

Recall from Appendix Al (see in particular Fig. 18) that the fundamental 
binary trees on (ra + 1) = 3 labeled leaves are of three types (twists are 
inessential). Choose for instance the tree T12 corresponding to the bracketing 
scheme 

T12 — ((JU2),,,J3)J (108) 

and pay attention to the fact that in the drawings we further simplify labels 
by setting ji = 1; J2 = 2; ji2 = 12; J3 = 3 while J is left unchanged. 
Since now we have always placed the labels onto the nodes (following the 
notations of (Topic 12)) but here we switch to the conventions of |T2] 
by labeling edges. Thus the tree acquires a new edge springing from the 
vertex formerly labeled by J and becomes a 3-valent graph -denoted by ti2- 
as shown in Fig. 6. We may interpret the vertices of such new graphs as 
"interaction vertices " for the pair of incoming spin variables. 

Figure 6 

The procedure to pass from ti2 to a decorated 2-disk with marked points is 
carried out in a few steps illustrated below. 

Step 1. A topological transformation can be performed on ti2 (tech- 
nically it corresponds to the action of a thickening functor, see e.g. 02]) 
which consists in "blowing inside" the graph and smoothing the corners. 
The resulting topological 2-manifold, depicted in Fig. 7, is a 2-sphere S'^ 
with boundaries represented by four circles (a "punctured" sphere). We la- 
bel such boundary circles by arbitrary labels {a, b, c, d} and denote such a 
surface by {S^; a, b, c, d). 

Figure 7 

Step 2. Fill up (5*^; a, 6, c, d) with an open set C M^: the resulting 3- 
space is (a portion of) a "precursor" of a handlebody M.^ (a handlebody 
is a closed oriented 3-manifold -the complement of a link in the 3-sphere 
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S^- obtained by gluing tubular neighbourhoods of the link along the bound- 
ary 2-disks). Here we simply glue closed 2-disks along the boundary circles 
of (5^; a, b, c, d) to get a 3-manifold with boundary -topologically equiva- 
lent to the 3-ball bounded by a 2-sphere- denoted by {D^, S*^) (a,b,c,d), where 
(a, b, c, d) may be thought of as labelings for 2-disks. The corresponding 
picture can be visualised by looking at the configuration in Fig. 7 as rep- 
resenting a "solid" object with 2-dimensional disks placed over the circles 
a, 6, c, d. 

Step 3. Consider an oriented embedding 

I, : ti2 P',5')(aAc,d) (109) 

which can be realized as a smooth map by splitting the 3-valent vertices in 
ti2 according to the rule illustrated in Fig. 8. 

Figure 8 

Then the image of ti2 into (-D^, S"^) (a,b,c,d) obtained by requiring that 1 ^ 
a, 2 ^ b, 3 ^ c, J ^ d as shown in Fig. 9 (we agree to drop out the 
auxiliary labels on the boundary 2-disks once the map it in (jl(J9|) has been 
implemented) . 

Figure 9 

Having started from the fundamental binary tree of p08|) we end up with 
a set of three (embedded) disjoint strands with crossings. In TQFTs con- 
figurations of this type represent "precursors" of knots/links, which are the 
observables of the theory as explained in the introductory remarks of this 
section. 

We may note that the procedure outlined so far depends on the choice of 
the embedding map ()109|) and thus it is not uniquely defined. As discussed 
in 1321 (Ch. 7) is always possible to establish a bijection between TQFTs be- 
longing to the categories FCx^^^ (graphs) and SDif f^^^ (smooth surfaces) 
but the construction is not well defined since one should actually consider 
"categories of all possible choices of embeddings and regular neighbourhoods, 
and get an induced TQFT by taking inverse limits over these categories". 
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Step 4. By taking "time slicings" on the configuration of embedded 
strands in Fig. 9 equipped witli tlie downward orientation we get tlie picture 
sliown in Fig. 10, wliere points and boundaries of disks inlierit consistent 
labelings. 

Figure 10 

Tlie final configuration (at the bottom of Fig. 10) represents a 2-disk (D^; 
ii! J25 jsi ji2', J) with three marked points and one marked circle inside (binary 
marked 2-disk for short). Note that 

• S'f/(2)-labelings of marked points (1,2,3) = {ji,j2,j3) can be freely 
chosen; 

• the labelings of the internal circle (12 = ji2) and of the boundary of the 
disk (J) have suitable ranges (c/r. (jH}) and are induced by the original 
tree structure. 

Step 5. If we take into account all the three fundamental binary coupling 
trees (each mapped into a suitable binary marked 2-disk) we realize that the 
algebraic structure relating the associated Hilbert spaces is encoded into 
the Racah identity (see the explicit expression (@7j)) which we write down 
schematically (apart from weights/phases) as 

'R- (Jl2, J3l) = (^'12' ^23) 7^ (J23, J31) (110) 

323 

The pictorial representation of piO|) acting on (the Hilbert spaces of) binary 
marked 2-disks is given by the triangular commutative diagram of Fig. 11 
which has the same content of the triangular graph shown in Fig. 22 of 
Appendix A2 (where the edges are thought of as topological moves on binary 
trees). 

Figure 11 

Step 6. By taking into account the different binary coupling schemes of 
{n + 1) =4 angular momenta we would get a combinatorial picture based 
on binary 2-disks with four marked points and two marked circles, each as- 
sociated with its own computational Hilbert space defined in (fT3j) . Then the 
natural algebraic structure linking the five spaces (up to phases) is provided 
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by the Biedenharn-Elliott identity (j46p encoded in the pentagonal diagram 
depicted in Fig. 23 of Appendix A2. 

Remark 6.2. As is well known from representation theory of simple Lie 
algebras (see e.g. [SI]), the multiple tensor products of irreducible modules 
can be handled formally by employing intertwiner spaces. More precisely, 
if A, A', A",... label highest weight representations of the algebra, the iso- 
morphisms Va ® Va./ = Va./ ® Va between the modules supporting (the tensor 
product of) the irreps A, A' is reflected into an isomorphism of intertwiner 
spaces 

F : T^^;, (111) 

and the three-fold isomorphism (Va®Va')®Va." = Va(8>(Va'®Va.") corresponds 
to an isomorphism 

i i 

of suitably defined intertwiner spaces. It can be shown that in fact the iso- 
morphisms ()111|) and pi2|) are all we need to treat arbitrary tensor products 
of (a finite number of) irreducible modules provided that three compatibil- 
ity conditions are fulfilled, namely a so-called pentagon and two hexagon 
identities. In the case of the Lie algebra sl{2) the pentagon relation is the 
Biedenharn-Elliott identity (j46p while the two hexagon relations become 
identical and coincide with the Racah identity ()47|) . Thus we recover the 
content of the Biedenharn-Louck Theorem of Section 4. 1 stated on the basis 
of SU (2) recoupling theory. A 

Summing up, the algebraic structure underlying the kinematics of the 
spin network simulator encodes automatically the pentagon and hexagon 
relations without resorting to ad hoc hypotheses. We have also shown that 
this combinatorial model for computation can be mapped -not uniquely- 
into the topological approach jlj and is not affected by locahsation problems 
typical of any "purely topological" setting. Finally, all the gates appearing 
in the spin network framework are unitary while unitary representations of 
the braid group must be carefully picked up to fit with the usual (Boolean) 
quantum circuit model (see |53,). For the convenience of the reader we collect 
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below a concise dictionary of the basic objects, spaces and maps employed 
in the two approaches. 

We argue that the partition functions of the two models may be related 
to each other much in the same way as the (regularized) Ponzano-Regge 
functional corresponds to a double Chern-Simons (c/r. Remark 5.1 at the 
end of Section 5 and in particular (j92|) ). 





combinatorial approach 


topological approach 


information 
encoded into 


binary coupling trees 
on 3 leaves 


topological qubits 

(^2 ^3pts) 


computational 
Hilbert spaces 




7^cs3(^',3pts) 


gates 


phase and Racah transforms 


unitary actions of i?3Ar 


compatibility 
conditions 


B-E (pentagon) identity 
Racah (hexagon) identity 


Yang-Baxter equation 



Table 1: A dictionary containing the basic ingredients of the combinatorial 
and topological approaches. 
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Appendix A. Spin network combinatorics 



The next two paragraphs are inspired by the basic reference jT^ (Topic 12) 
where binary couphngs of SU{2) angular momenta and unitary transfor- 
mations between pairs of such schemes (recouphng coefficients) are explored 
on the basis of their underlying graphs combinatorics. The origin of such an 
approach based on graph theory dates back to the Russian school of nuclear 
physics ([55j and earlier references therein) and gave rise to diagrammatical 
methods of vast applicability ^21; ^M- F^ck and collaborators have recently 
discussed some improvements concerning the efficiency of calculations for 3nj 
symbols [22] j [2S]- The latter achievements are summarised in Appendix A3, 
together with other results exploited in Section 4.3 in connection with spin 
network computational complexity. 

Al. Binary coupling trees 

A rooted binary coupling tree on (n + 1) leaves (terminal nodes) is a tree T 
-namely a connected graph with no cycles (closed loops)- characterised as 
follows. 

i) There exists a special vertex, the root. A 

ii) The tree is binary, namely each of its nodes has zero or two "siblings". 
More precisely, if we draw the tree with the root at the bottom and the leaves 
at the top as in all Figures of this section, the siblings of a particular node are 
the nodes lying in the nearest upper level which are connected to the given 
node by an edge. Thus the number of siblings (the out-degree) of both root 
and internal nodes is two, while each leaf has out-degree zero. The number 
of internal nodes of a rooted binary tree on (n + 1) leaves is {n — 1) and the 
total number of nodes is 2n+ root = 2n + 1 which of course coincides with 
the cardinality of the tree as a graph. A 

iii) The leaves are decorated with (n + 1) distinct labels. 

Generally speaking, we may attach to the leaves labelings {ii, 12,13, ■ ■ ■ , 
in+i} thought of as a permutation of the integers {1, 2, 3, ... , n+1}. However, 
in the framework of the quantum theory of angular momenta, labels are to 
be interpreted as quantum numbers {ji-^, ji^, ji^, . . . , ji^^-^} associated with 
a set of {n + 1) mutually commuting angular momentum operators Ji, J2, 
J3, . . . , J„+i. Once assigned any such a labeling to the leaves, we induce a 
consistent decoration on the other nodes by associating with them -moving 
downside along the tree- the quantum numbers of the intermediate angular 
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momentum operators arising from the pairwise couplings described by the 
tree. According to this rule the root acquires a label J , the quantum number 
of the total angular momentum J = Ji + J2+ ■ ■ ■ + J„+i. Referring to © and 
(|7j) of Section 2 it follows that rooted binary coupling trees on (n + 1) leaves 
represent in a faithful combinatorial way the structure of the computational 
Hilbert spaces ?i;J(b) given explicitly in (fT^ . 

As an example we sketch in Fig. 12 the rooted binary coupling tree 
corresponding to the particular binary bracketing structure of (fT^ or (jHJ.A 

Figure 12 

Remark A.l. As a matter of fact there exist other discrete structures in 
one-to-one correspondences with rooted binary trees and commonly used in 
theoretical computer science. Suppose you have a finite set of symbols S {e.g. 
letters in an alphabet A = {x,y, z, . . .}) endowed with a binary operation 
denoted by paired round parentheses 

{x,y) e (SxS) ^ (xy) (113) 

In the most general case the operation is neither commutative [{xy) 7^ {yx)] 
nor associative [{{xy)z) 7^ {x{yz))]. Note that we may enlarge the set 5" 
to include also and ")", but then we should require that the admissible 
words possess an equal number of left and right parentheses. As a particular 
example consider the case where S" = {(, )}, namely no other symbol except 
the parentheses themselves. The objects we obtain - with n and n ") " - 
are called Dyck words of length 2n [e.g. (()()) and ((())) for n = 3] and are 
enumerated by Catalan numbers. More generally, the structure induced by a 
binary operation on a finite set 5* turns out to be faithfully encoded in binary 
trees. As shown in Fig. 13, we may represent any pairing of two symbols {xy) 
by labeling with x and y, from left to right, two nodes in the tree which meet 
in a third node, then pairing the resulting symbol with a. z ^ S attached to 
another node, and so on (note that this procedure generates automatically a 
root in the tree). A 

Figure 13 

For the convenience of the reader, we collect here some useful facts from 
the combinatorics of rooted binary trees {erf. [13], |231 IHI]; IM] and 
references therein). Having in mind the definition of rooted binary tree given 
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in i)-ii) we are going to distinguish between "plane" and "not plane" unla- 
beled trees (here we use the adjective "plane" instead of "ordered" to avoid 
confusion with possible ordering of labels to be assigned to leaves: for the 
moment there is no labeling at all). A tree T is plane if its nodes -except 
the root- are put into an ordered partition of disjoint {Ti, T2, T3, . . . , T^} in 
the Euclidean plane, where each is a plane tree. In other words, in plane 
trees we distinguish between left /right and left /right subtrees. In Fig. 14 all 
plane unlabeled trees on 3 and 4 leaves are depicted. 



These trees (on (n + 1) leaves) are enumerated by Catalan numbers 



The first few terms in the sequence (jll4p are collected in Table 3 at the end 
of this paragraph, while more terms are listed in jHS] (ID Number A000108). 
According to Remark A.l, an identical counting holds true also for Dyck 
words and - as illustrated in Exercise 6.19 of (21] ^ there is actually a plenty 
of discrete structures from many branches of mathematics whose enumeration 
involves Catalan numbers. 

Not plane trees are obtained as equivalence classes of plane trees of the 
same size under reflections with respect to vertical axes through each node 
which is not a leaf. The resulting structures are known as "types" (of not 
plane trees) and enumerated by Wedderburn-Etherington (W-E) numbers 
Bn+i for which a closed form expression is not known. The first terms of 
this sequence are listed in Table 3 (see jSH]; ID Number A001190 for more 
terms). In Fig. 15 the trees of W-E-types on 3 and 4 leaves are shown. 



The classification for unlabeled trees considered so far is summarised in the 
upper row of Table 2 at the end of the section. 

We say that a tree is labeled when we assign some symbol to each of its 
leaves. Generally speaking, we may use either distinct labels (as in binary 
coupling trees, see iii) above) or a binary label {0, 1} (as in search trees) or 
even a same label for each leaf. 



Figure 14 




(114) 



Figure 15 
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Remark A. 2. To illustrate the subtleties arising from assignments of no 
label, a same label or distinct labels to leaves of trees, let us consider Catalan 
trees again. According to our previous definition they should be unlabeled, 
but it is easily recognised that the same combinatorics is shared by two more 
plane binary tree structures (see Fig. 16), namely 

• Trees with all leaves labeled by a same x. In this case the trees encode 
a non associative, partially commuting binary operation on the alpha- 
bet {x}. ["Partially commuting" means that inside each pairing {xx) 
commutativity is trivially ensured, but for instance {{x^)x) 7^ (x(a;^)]. 
Note that the case of a non associative totally commuting binary oper- 
ation corresponds to W-E-types labeled with a same x since we would 
get in that case = as should be clear by looking at Fig. 
15, top). 

• Trees with a fixed sequence of distinct labels {e.g. ordered lexicograph- 
ically from left to right) on their leaves. 

In the following we shall refer to the above labeled structures as Catalan 
trees too. A 



Figure 16 



From now on we are dealing with distinct labelings, chosen for simplicity 
in a string of (n -|- 1) Latin letters {a, b, c,d, . . .}. 

Starting with the plane category, we have to decorate the leaves of Catalan 
trees with all possible permutations of labels. Thus the number of these 
objects for any n is given by (n + 1)! C„, where C„ is the Catalan number 
defined in flll4|) . This number, usually written as 

C. ^ <^ . (115) 

nl 

is the quadruple factorial and in Fig. 17 all labeled plane trees on n -|- 
1 = 3 leaves are shown {C2 = 3!C2 = 12). The first terms of this integer 
sequence are listed in Table 3 and more terms can be found in [2^1; ID 
Number A001813. 



Figure 17 
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To enumerate labeled not plane trees on (n + 1) leaves we realize that 
there exist n axes through the (n — 1) internal nodes plus root. Thus, on the 
basis of (jllSp and looking also at Fig. 17, we have to drop out exactly 2"^ 
configurations since the trees are binary (namely each of the former n nodes 
has exactly two siblings). The resulting counting reads 

where (2n — 1)!! = 1 ■ 3 ■ 5 ■ 7 . . . is the double factorial (see Table 3 and 
ID Number A001147). The combinatorics of this enumeration can be also 
understood by picking up each W-E type of unlabeled trees (see Fig. 15) 
and decorating consistently its leaves. Then 

b=l 

where Pn+i (b) represents the number of ways to decorate leaves of type b 
trees with (n + 1) distinct labels in such a way that (xy) = (yx) for each 
binary parenthesization involving both nodes and subtrees. In Fig. 18 not 
plane labeled trees on 3 leaves are depicted {D2 = C'il'^ = 3): they arise 
from the unique W-E type shown in Fig. 15, top. 

Figure 18 

In the following Table 2 the combinatorial enumerations for all trees con- 
sidered so far are summarised, while in Table 3 the first few terms of the four 
integer sequences are written down to give the reader an idea of their rates 
of growth. 
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plane 


not plane 


unlabeled 


Catalan number 

(2n)! 
^" n\(n+l)\ 


Wedderburn-Etherington number 

Bn+l 


labeled 


Quadruple factorial 
r - (2'^)! 


Double factorial 
D„ = g = (2n-1)!! 




= card {0n{V, E) 


= card(6„(V,E) 



Table 2: Enumerations of rooted binary trees on (n + 1) leaves according 
to plane/not plane and unlabeled/labeled categories. In case of labeled 
trees the numbers represent the cardinalities of the Twist-Rotation and Ro- 
tation graphs, respectively (see Appendix A2). 



n 


1 


2 


3 


4 


5 


6 


7 


Bn+l 


1 


1 


2 


3 


6 


11 


23 


Cn 


1 


2 


5 


14 


42 


132 


429 


Dn 


1 


3 


15 


105 


945 


10395 


135135 


Cn 


2 


12 


120 


1680 


30240 


665280 


17297280 



Table 3: The first terms of the four integer sequences of Table 1 arranged 
according to their rate of grow. -Bn+i are Wedderburn-Etherington numbers; 
Cn are Catalan numbers; Dn are double factorial numbers; C*„ are quadruple 
factorial numbers 

A2. Twist-Rotation and Rotation graphs 

Here we explain -keeping on using extensively graph-theoretical tools as in 
[ISI, 122]^ the construction underlying the discrete computational space of 
the spin network simulator discussed in Section 4.1. For simplicity we do not 
change the notation used there, although many definitions and results can 
be applied also to larger classes of graphs (the so-called "distance graphs" 
see ISZl)- 

Given a rooted binary tree T on (n + 1) labeled leaves {n > 2), two 
kinds of topological operations (moves) can be considered, namely rotations 
and twists ("rotation" is not to be confused with rotation matrices or op- 
erators used in Sections 3.2, 4.1 and Appendices B). These moves represent 
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alterations in the shape of the tree generated around either an internal node 
(rotation) or a node which is not a leaf (twist) and all trees in a given class 
can be reached on applying admissible moves to an arbitrary tree chosen in 
that class (roughly speaking the moves are "ergodic"). 

Rotations around non-root internal nodes consist in swapping subtrees 
or nodes as shown in Fig. 19. In the upper region there appears on the left 
a generic (portion of a) tree T with subtrees A, B, C and a fourth subtree 
R containing the root: the rotation around the node x transforms T into T' 
and, conversely, the rotation around x' G T' changes back T' into T . In 
the lower region a rotation on a particular tree is shown, together with its 
inverse operation. 

Figure 19 

Since there are {n — 1) internal nodes, {n — 1) different rotations can be per- 
formed on any tree (with fixed labelings on its leaves). Note that a rotation 
alterates the shape of the tree near the corresponding node but leaves the 
rest of its structure intact. 

Twists around non terminal nodes consist in exchanging left and right 
subtrees (or nodes). We draw in Fig. 20 both a general twist around a node 
X & T and a particular twist on a specific tree. 

Figure 20 

There are n possible twists on a rooted labeled tree on (n + 1) leaves and 
these transformations may alterate the global shape of the tree. 

Coming back to combinatorics of labeled rooted binary trees (see Table 
2), we realize that every category may support rotations but only "plane" 
trees (enumerated by the quadruple factorial Cn) admit twists. For what 
concerns labeled Catalan trees in Fig. 16 (bottom) we see that they are 
actually connected by the rotation represented in Fig. 19, but no twist is 
allowed since the sequence of labels must be kept fixed. Looking at Fig. 
17, where plane trees labeled in all possible ways are shown, we see clearly 
that both rotations and twists appear as admissible moves connecting pairs. 
Finally, not plane trees (enumerated by the double factorial Dn and shown 
in Fig. 18) undergo only rotations since they do not distinguish between left 
and right at each non terminal node. The further step consists in building 
up new graphs associated with each of the above categories (these structures 
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are called "distance graphs" and the reason for such terminology will become 
clear in the next paragraph, where distances will be introduced). Denote 
generally such a graph by G„(V, E) (as in Section 4.1), where V and E are 
the vertex and the edge sets, respectively. The vertices are in one-to-one 
correspondence with rooted binary trees on (n + 1) labeled leaves, namely 

t(") < — > V eV, (118) 

where T*^""^ stands for a tree in a given class. The edge set E is generated by 
linking pairs of vertices v,v' G by an undirected edge e{v,v') if, and only 
if, the corresponding trees are related to each other either by a rotation or 
by a twist. Formally 

e{v,v') eE^T^-^ ^ T^;\ (119) 

where ^ stands for one of the topological moves described above. According 
to the remarks on admissible moves acting in different categories of trees we 
speciahse G„(V,i?) as follows. 

• Labeled plane trees on {n + 1) leaves represent the vertex set of the 
Twist-Rotation graph Gj^^{V,E) denoted as in (jSZj), namely 

Gl^{V,E) = en{V,E). 
Its cardinality is the quadruple factorial number given in ()115|) 

card(0„(\/,E)) = \V\ = Cn (120) 

and each edge represents either a rotation or a twist (c/r. also ()38|) ). 

For any n > 2 &n{V, E) is an undirected regular cubic graph (namely 
there is no ordering on its vertices and each vertex has valence three). It 
is also a planar graph, i.e. it can be drawn onto the 2-dimensional sphere 
without crossings. For n = 2 card(C52) = 12 (c/r. the trees in Fig. 17) 
and 02 has the shape of a truncated tetrahedron made up of triangular and 
hexagonal faces. 03 has cardinality 120 and represents the graph of a 3- 
valent polyhedron made up of pentagons and hexagons. A portion of this 
graph is shown in Fig. 21 (reprinted from j22I) and its shape looks quite 
familiar since the discover of fuUerene. 

Figure 21 
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• Labeled not plane trees on (n + 1) leaves represent the vertex set of 
the Rotation graph G^{V,E) denoted as in (jSHI), namely 

G^{V,E) = 0n{V,E). 
Its cardinality is the quadruple factorial number given in ()116|) 

card(6„(y,E)) = IV^I = Dn (121) 

and each edge represents a rotation (c/r. also (jlOI))- For any n > 2, &niV, E) 
is an undirected regular graph of valence 2(n — 1) which turns out to be not 
planar for n > 3. 

The graph (&2 with cardinality 3 shown in Fig. 22 is a trivial example 
from the combinatorial point of view but - by exploiting the encoding map 
(|78p - we see that it shares the same content as the "triangular" Racah 
algebraic identity (jTTj) introduced in Section 4.1. 

Figure 22 

The graph 03 (whose vertices are trees on 4 leaves) is depicted in Fig. 1 
of Section 4.1: from a topological point of view we see that some vertices 
have been doubled to avoid crossings; when we identify by an antipodal 
mapping opposite vertices and edges we realize that this graph lies in the 
real projective space MP^ (namely it is not planar). Other pictures of 03 can 
be found in [22] and |2S1- In Fig- 23 one of the pentagons belonging to both 
03 in Fig. 1 and 03 in Fig. 21 is depicted: edges represent rotations from 
a topological point of view and at the same time the figure is the diagram 
encoding through (fTHj) the Biedenharn-Elliot (pentagon) identity (jiUj) . 

Figure 23 

• Labeled Catalan trees on [n + 1) leaves represent the vertex set of 
the Rotation graph denoted by 

G^{V,E) = Q^{V,E) (122) 

to distinguish it from the previously defined 0n. QniV, E) is a regular graph 
of valence (n — 1) and its cardinality is given by the Catalan number ()114|) . 
namely 

card(0„(\/,E)) = \V\ = C„. (123) 
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Although this kind of Rotation graph is not suitable to model the com- 
putational space of the quantum simulator (Catalan coupling trees are not 
sufficiently general) we meet these structures in connection with issues on 
combinatorial complexity discussed in Section 4.3 and in the next paragraph. 

A3. Combinatorial complexity 

All graphs described in Section A2 are actually "distance graphs" |57j since 
they encode in their edge sets operations on the basic objects (trees) as- 
sociated with their vertex sets. To evaluate quantitatively how far away 
pairs of rooted labeled binary trees are -and using standard terminology in 
discrete mathematics- we introduce explicitly Twist-Rotation and Rotation 
distances. 

Given two binary trees T^^\ t!^^ -thought of as vertices in (S„, 0„, g„, 
respectively- their distance is the length of the shortest path joining them, 
namely the minimum number of topological operations needed to transform 
one tree into the other. With an obvious meaning of symbols we set 

min {lengths of paths T^'^ Tg"''' in 0„} 
min {lengths of paths ^ Tg"'' in 
min {lengths of pathsT^*'"^ ^ t!^'^'' ing„}. (124) 

The diameters of the above graphs are naturally defined in terms of distances 
according to 

Dzam(60 = max {d™(T^"\Ti")) I T("\Tf) G 6 J 
Z^mm(0„) = max{rf^(T("\Ti"^)|rf"\T(") G 
Diam (0„) = max {<r (T^"), T^^) | T^"), T^^ G gj. (125) 

In addressing combinatorial complexity we are primarily interested in com- 
puting or estimating such distance functions and diameters since in Section 
4.3 we relate combinatorial complexity to computational complexity of the 
spin network simulator on the basis of the quantum encoding maps intro- 
duced in fITBj) and (fTSj) . In this respect it may be useful to collect here some 
known results from graph theory (addressed in a classical computational 
complexity framework). 



ci^(Tf"\Ti")) = 
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• For what concerns Catalan trees equipped with rotation distance d'^, 
a major breakthrough was achieved in jHH], where the authors proved the 
existence of a tight bound on the diameter of g„ given by 

Diam{gn) < 2n - 8 (126) 

for trees on (n + 1) terminal nodes. Their elegant proof relies on the trans- 
lation of the combinatorial problem for trees into an equivalent geometrical 
one, namely enumerating triangulations of a polygon with {n + 1) edges and 
finding the maximum number of diagonal flips needed to convert one trian- 
gulation into another (diagonal flips are topological moves on 2-dimensional 
triangulations which transform a quadrilateral dissected into two triangles 
into the configuration generated by cutting along the other diagonal). The 
argument to get the final result involves the construction of 3-dimensional 
hyperbolic polyhedra and calls into play volume estimates in hyperbolic ge- 
ometry. 

• In [23, purely combinatorial counting tools used previously in ^S] are 
improved to establish lower and upper bounds for the diameter of the Rota- 
tion graph (3n- The explicit form of the upper bound is (Ig = log2) 

Diam{(5n) < n\g{n) + n ~ 2\g{n) + 1 (127) 

which is compatible with older estimates (see e.g. |E3)- Therefore we may 
conclude that the diameter of C5„ grows no faster than polynomially in n, 
namely 

Diam{(5n) < n\g{n) + 0{n). (128) 

• A crucial open question remains the complexity status of computing the 
rotation distance even for the simplest case of on the graph g„ made up by 
unlabeled Catalan trees (binary search trees in computer science language). 
In particular (c/r. the recent papers IHI], [02]) it is not known 

• whether the problem is NP-complete; 

• whether d^ can be determined in time polynomial in n (that is to say, 
whether there exists an efficient algorithm to compute it exactly). 

There are however (classical) polynomial time algorithms which estimate this 
distance (or its lower/upper bounds) ([EI], IHSI)- 
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• The two remarks above hold for the rotation distance too, since 0„ 
is much bigger than g„: more precisely, if we think about "labeled" Catalan 
trees (as pointed out in Remark A. 2, Appendix Al) then 0„ is actually a 
subgraph of 0„. To our knowledge there is no algorithm for computing 

in (J5„, although we may consider lower and upper bounds on Diam [23] 
(which are both of order ?T,lg(n)) as an estimate of the number of elementary 
operations -Racah transforms- appearing in an "optimal" expansion of a 
3nj symbol. (Note however that such optimal expansion does exist only if 
we could actually computing d^{T^\T^'') for arbitrary T^\T!f''^ G S5ri) 

• In order to try to overcome the difficulties outlined above, a basic 
strategy emerges in current literature, namely the idea of introducing in 
Qn some more specific notion of distance function. In 64^ a "restricted" 
rotation distance has been considered, and linear lower and upper bounds 
are estimated in terms of the number of interior nodes of trees. In |^ 
a "right-arm" rotation distance is shown to be computable by an efficient 
algorithm (in polynomial time O(n^) for trees on {n + 1) leaves). 

• In the field of molecular biology, phylogenetic trees (or dendograms) 
for groups of species turn out to be powerful tools to address the study of 
similarities and dissimilarities appearing in biological evolution (see e.g. jHS] 
and references therein). Theoretical biologists use different kinds of defini- 
tions with respect to our previous classifications from discrete mathematics. 
In particular, their rooted binary phylogenetic trees may have labels only on 
leaves (the different species) while the role of internal nodes is secondary. 
The distance function which is commonly used in this context is the "nearest 
neighbour interchange" (nni) distance which is surely more restrictive than 
the rotation distances introduced in ()124|) . although it is not so easy to es- 
tablish connections with the restricted rotation distances considered in our 
previous remark. The authors of [UH] claim to have proved the fact that 
computing the nni distance is NP-complete (both in the unlabeled and in 
the labeled cases) by a reduction from Exact Cover by 3-sets) (X3C) which 
is known to be an NP-complete problem |47j . 
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Appendix B 



Bl. Composition of W— rotation matrices 

In this Appendix, following the standard reference we illustrate the 
derivation of the symbolic expression given in (j^^ which represents -in the 
language of quantum theory of angular momenta- the generalised Clebsch- 
Gordan expansion involving a finite number of W-rotation matrices with the 
same arguments. 

Denote as usual by 

Ji5 J2 ) J3 5 • • • 5 Jat; 
Ji,j2, is, ■■■Jn, 

mi,m2, ms, ...,mN, 

m\,m'^, m'g, ...,m'^ (129) 

an set of commuting angular momentum operators and the corresponding 
sets of quantum numbers. 
Consider the operators 

= Jl + J2 + J3 + ••■ + Ji (130) 

which are defined, for each 1 < i < A^, by taking any kind of binary coupling 
consistent with this vector addition rules. We use here /C's instead of K's 
used in Section 2 to stress the fact that we can choose anyone of the schemes 
considered there (note that in (jl3Up there appear also /Ci = Ji and /Cat = 
J, so that the counting of intermediate operators is not in contradiction with 
previous statements). The quantum numbers of these operators are denoted 
by 

Kl , /t2 , , . . . Ktv, (131) 
while (two possible sets of) magnetic quantum numbers are given by 

Mi= mi + m2 + ms + h mj 

M[ = m[ + m'^ + m's + ■■■+ m[. (132) 
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Coming to the composition of W-matrices, we start by considering the 
Clebsch-Gordan series, namely the expansion of the product of two W- 
matrices (with the same arguments) labeled by ji and j2- It reads 

D^' , (a/57) , (a/37) = 

^/imij2m2^MM'i'^Pl)Cfim{j2m'2^ (133) 

J=\ji-j2\ MM' 

where we denote simply by J the sum Ji + J2 and by J, M, M' the cor- 
responding quantum numbers (leaving aside for the moment the notations 
introduced in p3()|l . ()131|1 and ()132j) ). By using the orthogonality condition 
of the Clebsch-Gordan coefficients C" „ we can invert p33p to get 



Dmm' («/?7) = 

E C/Znm2 D'mim', M7) D^^^^, (a/37) C^i^^ ,,m^ , (134) 



mi m2 



where the j's entries of each Clebsch-Gordan coefficient must fulfill the tri- 
angular inequality (here in particular we have just to require \ji — < <^ < 

Ji + J2)- 

The decomposition established above represents the starting point to ad- 
dress the more general case of > 2 incoming angular momenta. By fully 
restoring the notations explained at the beginning of this section and apply- 
ing successively (|134p together with the orthogonality conditions whenever 
necessary, we get the generalised C-G expansion involving the product of N 
W-rotation matrices. Its explicit expression reads 



E fl C::^k..nm. D-^m^ i^h) C:%_^ , (135) 



mi...m]v 1=1 
m[...m'j^ 



where it has been assumed kq= Mq = Mq = and all the triads {ji , , Ki} 
satisfy triangular inequalities. 
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We may recast the previous expansion in a symbolic form ( cfr. (I29|) in Section 
3.2) by keeping only the dimensions of representations both in D's and C's 
(the typographical changes are introduced to remind that these objects are 
matrices) 

N 

E n {C^l.. ^"M7) C^t,,). (136) 

{m, m'} i=l 

Here we set 

J = Kn = Ji + J2 + J3 + 1- Jat; 

J =kn (137) 

to comply with the notation used in the main text and in the particular 
case ()133|) (consequently the matrix indices oi (0/^7) are intended to be 
M = mi + 772.2 + 'TT'S + ■ ■ ■ + w^at; M' = m[ + m'2 + m'g + ■ ■ ■ + m'^). 
The summation over the magnetic numbers of the incoming J's which ap- 
pears in front of ()135j) (or ()136j) ) is a partial trace on the iV-product of 
triples contained in round parenthesis: any such triple combination can be 
directly evaluated by substituting the numerical values of the "elementary" 
W-matrices and of the suitable pair of C-G coefficients. 

The structure displayed in ()135p is drastically simplified if we look at two 
particular cases. 

• Fermionic case (A^ spins = | in the symmetric multiplet) 
For 

3i = 32 = k = ■ ■ ■ = 3n = 1/2 

and 

Ki+i = Ki + 1/2 J = kn = N/2) 
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the expansion (|135p becomes 



N 



A -Di^, {a(3^) = E n ^Im' («/57) (138) 



r?ii+...+mjv=M 4=1 
/ _ 

JV" 



m[+...+m'j^=M' 



where the weight A is given by 

v/(J + M)! (J - M)! (J + M')!(J - M')\ / (2J)! (139) 

and each factor -D^m' represents the W-matrix in the fundamental represen- 
tation written in terms of the Euler angles a /3 7, namely 



^mm' l«P7J - ^ gia/2gi^(^/2) e-*^/2 e*«/2 cos(/3/2) 6^^/^ j- U^uj 



• Bosonic case [N spins = 1 in the symmetric multiplet) 
For 

ji = j2 = js = ■ ■ ■ = Jn = l 

and 

Ki+i = + 1 J = Kiv = A^) 

the expansion ()135p becomes 

A ■D^.jvf M7) = 



TV 

5^ n V^(l + 5™.o) (1 + 5^0) («/?7), (141) 

mi + ...+m]v=M i=l 
m'^+...+m'^=M' 

where A is given in (jl39|) and each factor Dl^^, represents the W-matrix in 
the j = 1 representation written in terms of the Euler angles a jSj, namely 

D'mm' («/?7) = e-^*^" iP) e-""'' (142) 
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with 



/ l+cos/3 
' 2 



sin /3 



^ COS/3 
sin /3 



V2 

1— cos /3 

2 



l-cos/3 \ 
2 » 
sin /3 

1+cos /3 
2 



(143) 



Remark B.l. Equations ()138|1 and ()14ip turn out to be useful for evalu- 
ating in practice matrix elements of W-rotations acting on irreducible repre- 
sentation spaces. Notice in particular that matrix elements flMlj) can be actu- 
ally obtained by considering the symmetric combination of two spins = 1/2, 
namely by exploiting p38|) with N = 2. More generally, any (2jj + l)x 
{2ji + 1) W-matrix (as those appearing in ()135|l ) acts on a representa- 
tion space which is irreducible under the action of SU(2) and (according 
to Majorana splitting jHSl) such eigenspace may be always viewed as a com- 
pletely symmetric multiplet generated by 2jj | spin kinematically indepen- 
dent "particles". By taking advantage of this remark we may factorize each 
^, in terms of p4()j) . as discussed in relation with the estimate given in 



B2. t/— rotation matrices 

In this appendix we collect from some useful definitions and examples 
concerning rotation matrices expressed in terms of the direction of the rota- 
tion axis and of the rotation angles, namely alternative forms for the M-gates 
introduced in Section 3.2. 

Let n(0, $) a unit vector in the direction of the rotation axis and u the 
rotation angle. An element of an [/-matrix is defined according to 

Ui^M, {to; e, $) = {JM I exp{-zcjn ■ J} | JM' ) (144) 

in a representation space labeled by the quantum numbers J, M. The W- 
rotation matrices introduced in are related to ()144|) by 



[/j^^, (u;; e, $) = ^Di,^„($,e,-<l>)e-*^"-Di,„^,($,-e,-<l>). (145) 

M" 

For rotations around coordinate axes we find the expressions (the reduced 
VT-matrices df^j^^j, were introduced in and given explicitly in ()14U|) and 
dug) of Appendix Bl for J = i, 1). 
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(x axis; 9 = |; $ = 0): 



f/MM'(^;2'0) = ^MM'(2'^'-2) = (-^) 4/M'M (146) 
.(yaxis;e = f;$ = f): 

UmnA^;\.\) = ^MAf'(o,^,o) = 4/M'H (147) 

• {z axis; 6 = 0): 

Ui,j^,{uj-^,^) = D](,A,,(0,a;,0) = W'e-**'^" (148) 
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Figure 1: The rotation graph ^3 where each vertex represents a binary cou- 
phng scheme of (ra + 1) = 4 angular momenta. Vertices on the perimeter 
have to be identified through the antipodal map, showing that 03 is not 
planar. Using the notations of Appendix Al, the 1 : 1 correspondences be- 
tween vertices and binary parenthesization on four letters, say {a,b,c,d}, 
read: 1 ^ (rf(6(ac))); 2 ^ (6(rf(ca))); 3 ^ ((ac)(M)); 4 ^ (d(a(6c))); 5 ^ 
{d{c{ab))); 6 ^ (c(rf(a6))); 7 ^ ((a6)(cci)); 8 ^ (a(d(6c))); 9 ^ ((arf)(6c)); 
10 ^ (a(6(cc/))); 11 ^ (a(c(M))); 12 ^ (c(a(M))); 13 ^ (c(6(arf))); 
14 ^ (6(c(da))); 15 ^ (6(a(cd))). 
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Figure 2: The computational setting based on the combinatorics of the Ro- 
tation graph 63 can be pictorially represented by blowing out each 
vertex into the corresponding computational Hilbert space and by looking 
at the gates as elementary operations on the base space and/or on the fiber 
(c/r. the remark at the end of this section and the analysis carried out in the 
next section). 
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Figure 3: Left: a marked 2-disk corresponding to the Hilbert space 'Hcsz 
(Z^^,3pts) with {a, 6, c, d} G {0,1,2,3}. Right: the 2-disk supporting the 
topological qubit y(D^,3pts). 
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Figure 5: The Yang-Baxter equation -R12-R13-R23 = -^23-^13-^12- 
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Figure 7: The punctured 2-sphere (S*^; a, b, c, d) with labels on the boundary 
circles. 




Figure 8: An interacting 3-valent vertex with two incoming edges is splitted 
by requiring that the first entry in the binary coupling overcrosses the second 
entry (and viceversa). 
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Fi gure 9: An embedding of the fundamental graph ti2 into {D'^, S"^) (^a,b,c,d) 
where we drop out the labels (a, b, c, d) on the boundary circles which ap- 
peared in Fig. 7. 
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Figure 12: A particular rooted binary coupling tree on {n + 1) leaves where 
incoming angular momenta are coupled sequentially. 
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X y z w 



(((xy)z)w) 




X y z w 



((xy)(zw)) 




Figure 13: Explicit bijections between binary parenthesizations over an al- 
phabet of 4 letters and rooted binary trees on 4 leaves. Here the ordering 
of {x,y,z,w} is fixed and only two correspondences are depicted [the other 
ones arise from {x{y{zw))); {{x{yz))w); {x{{yz)w))]. 
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Figure 16: Presentation of Catalan trees on 3 leaves as labeled trees: a same 
label assigned to each leave (a); a fixed sequence of distinct labels assigned 
to leaves (b). 
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Figure 18: Not plane rooted binary trees on 3 leaves labeled with no distinc- 
tion between left and right. 
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a b c a b c 




Figure 19: A general rotation (top) and a particular rotation (bottom) on 
rooted binary trees. A,B,C are generic subtrees, and R is a subtree containing 
the root. 
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Figure 21: A portion of the Twist-Rotation graph ©3 representing the 
vertex-edge set of a truncated icosahedron. Only 30 of the 60 vertices are 
shown, and the picture can be completed by taking the mirror image of each 
tree at the antipodal vertex. The remaining 60 vertices of ^3 are arranged 
into an isomorphic graph which is obtained from the former one by exchan- 
ging a pair of labels, e.g. (a, b) — > (6, a). 
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Figure 22: The Rotation graph 02 is depicted as a triangle with binary 
trees attached to its vertices. Edges correspond to rotations (topologically) 
but - when deahng with binary couphng trees - they represent three Racah 
transforms related by the Racah identity (see also Section 4.1, Section 6 and 
Fig. 11). Similar configurations are associated to triangular plaquettes in the 
Rotation graph (S3 of Fig.l. 



21 



abed 




Figure 23: A pentagonal plaquette in the Twist-Rotation graph 03 of Fig.21 
or in the Rotation graph (Bs of Fig.l. When binary coupling trees are conside- 
red, edges represent Racah transforms related by the pentagon (Biedenharn- 
EUiott) identity, see Section 4.1. 
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